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Abstract. We study the phenomena of energy concentration for the critical 
0(3) sigma model, also known as the wave map flow from R 2+1 Minkowski 
space into the sphere S 2 . We establish rigorously and constructively existence 
of a set of smooth initial data resulting in a dynamic finite time formation 
of singularities. The construction and analysis is done in the context of the 
k-equivariant symmetry reduction, and we restrict to maps with homotopy 
class k 4. The concentration mechanism we uncover is essentially due to a 
resonant self-focusing (shrinking) of a corresponding harmonic map. We show 
that the phenomenon is generic (e.g. in certain Sobolev spaces) in that it 
persists under small perturbations of initial data, while the resulting blowup 
is bounded by a log-modified self-similar asymptotic. 



1. Introduction 

One of the simplest non-trivial models of Quantum Field Theory is based on 
the (2 + 1) dimensional Lorentz invariant 0(3) classical a-model. It is a nonlinear 
scalar field Lagrangian theory for a map <£> : R 2+1 — > § 2 C K 3 with the Lagrangian 
density: 

(1) £[«] = \d a <5> ■ m a0 , 

where m a p is the Minkowski metric on K. 2+1 . Evolution of the nonlinear scalar field 
<!> is described by the Euler-Lagrange equations: 

(2) □$ = -$(<9 Q $ • d a §) . 



The equation © belongs to the more general class of "wave-map" problems, in 
which <I> is a map from Minkowski space M 2+1 to a Riemannian manifold (M,g). 
The map $ is a solution of the Euler-Lagrange equations: 

(3) D a d a <S> = , 
corresponding to the Lagrangian density: 

(4) £[$] = X - 9l , d a &8p& m af} . 

Here {3> 1 } denote local coordinates on A4, which in turn (under the map) depend 
on the Minkowski variables {x a } a =o,i,2- D is the pullback of the Levi-Civita con- 
nection to the (trivial) bundle Q*(TA4). In terms of the local coordinates this 

This work was conducted when I.R. was visiting the Department of Mathematics, MIT. He 
was also partially supported by the NSF grant DMS-0406627. J.S. was supported by an NSF 
postdoctoral fellowship. 
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pull-back connection acting on sections of $*(TA4) reads: 

(5) D a = d a +T k aj , T k aj = T^)8 a ^ , 

where 1^ is the Christoffel symbol in the coordinates The wave-map equation 

§3§ then has the intrinsic form: 

(6) d a d a <$> k = -T k ;A<5>)d a <$> l d a & . 



The goal of this paper is to establish and rigorously analyze a catastrophic in- 
stability in the (2 + 1) dimensional 0(3) a- model represented by the equation ((2|). 
We will exhibit a spontaneous and monotonic self-focusing mechanism responsible 
for a dynamic formation of singularities for a rather large and stable set of initial 
data. This will be done through an entirely explicit and constructiv^ description 
of this phenomena. Our basic result is as follows: 



Theorem 1.1. For every < e <C 1 and 4 ^ k there exists a set of smooth initial 
data ($0)^o) G (§ 2 ,T§ 2 ) with energy E = 4irk + e 2 (i.e. the Dirichlet energy 
defined below), and a finite time T** = T**(<i>o, $0), such that the corresponding 
solution $(t,x) of problem @ remains smooth on the interval [0, T**) and develops 
a singularity at T** . More specifically, there exists a (smooth) decomposition $ = 
$ + TZ, such that as t — > T** we have that for any large < M an L°° bound of 
the forn^: 



M ^ vIMr** - 1)\ 

as well as a uniform bound on the energy of the remainder: 
(8) E[R] < e 2 . 

Furthermore, sufficiently small equivariant perturbations of ($o,$o) a ^o lead to 
blowup with the bounds ([7]) -(|5]) . 



The problem of a finite time breakdown of solutions of the problem @ has been 
a subject of intense study. From a purely analytical perspective, the context is the 
global regularity theory for the general wave-map equations (J3]), where it is sus- 
pected that the formulation of singularities is ultimately tied to certain convexity 
properties of the target manifold M. . 

From a more physical or gauge theoretic perspective, and in a specific context 
of the 0(3) model, the issue of possible singularity development is thought to be 
connected to the incompleteness of a certain moduli space which characterizes the 



This is in contrast to some of the examples of the focussing nonlinear Schrodinger and wave 
equations, where a finite time blow-up can be shown by non-constructive arguments (see 1131 and 
[25]). See however the work of Martel/Merle 1271 on the critical KdV problem for an example of 
a constructive finite time blow up mechanism. 

2 More precise asymptotic behavior in terms of the energy concentration will be given below, 
including both upper and lower bounds. 
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associated static solutions, and provides an approximation for the dynamical evo- 
lution through the so called geodesic hypothesis. Due to its analytic and physical 
interest, the equation @ has also been a popular subject of numerical and heuristic 
studies which all universally pointed in the direction of singularity formation. 

Before proceeding with a more detailed description of our main result, we believe 
it is useful to give a more thorough description of these various points of view. This 
begins with a discussion of the static solutions of @, that is harmonic maps into 
the sphere. Historically, one of the primary motivating factors of interest in the 
0(3) cr-model was due to the richness of the set of its static solutions. An ingenious 
procedure of Belavin and Polyakov [1] allows one to find these solutions in a given 
homotopy class characterized by the topological degree: 



k = — I $*(dA&) , 
4tt J r2 



as solutions of the first order "Bogomol'nyi equations" (see [4]). To realize this, 
one factors the energy functional: 

(9) V[$] = ~ f V K $ • V x $ dx , 
which is the potential part of the Dirichlet type energy: 

(10) £[$](*) = \ [ (d t <P-d t <S> + V x <P-V x $) dx = T[$\(t) + V[$\(t) , 
associated with the action of the Lagrangian |T]) : 

(11) / dxdt = - { (T[$] -V [$])(*) dt . 

Using the notation for the antisymmetric tensor on two indices, this factorization 
reads: 



(12) 



(9,$ ± e i J $ x 9,-$) • (9*$ ± x 9,-$) 



dx 



± \ ( ' {di$ x dj®) dx 

2 Jr2 



{d^ ± e- J $ x • ± e ij $ x 5,$) 



(ice ± 4nk 



from which it is more or less immediate that an absolute minimum of the energy 
functional V[$] in a given topological sector k must be a solution of the equation: 

(13) 9 i $±e i i $ x 9,-$ = . 

In terms of complex coordinates on R 2 and S 2 , the identities p3|) are seen to be 
nothing other than the Cauchy-Riemann and conjugate Cauchy-Riemann equations 
(this is a general phenomena, see [S]). Therefore, the moduli space of static 
energy minimizing solutions of |(5J| in a homotopy sector k can be identified with 
the rational maps (in z or z using complex variables) / : C — > C with degree k. Of 
particular importance to us will be the k-equivariant static solutions (of positive 



polarity) which are defined via the relation I(e l6 z) 
solutions by I k and refer to them as solitons. 



J6k 



I(z). We will label such 
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Having uncovered the structure of the space of static minimizing solutions, one 
is then led to the "geodesic" ansatz alluded to above for the approximate dynamics 
of time-dependent solutions (see also [25]). To understand this, the first thing 
to notice is that minimizers of the variational problem ([9]) are highly degenerate. 
Specifically, they are invariant under the full conformal group of linear fractional 
transformations acting on C. If one restricts to £;-equivariant solutions, then most 
of this symmetry is lost, and the only remaining degree of freedom which fixes the 
energy © is the scaling transformations: 

(14) I k {t,x) -> I k (t,x) = I k (Xt,Xx) . 

Based on this, one would expect that the path of least resistance according to the 
action (jlip would be for (symmetric) solutions close to the family of static solutions 
to "slide" along the moduli space Mk via the transformations (fT"4")) . That is, for a 
fully dynamic solution $ sufficiently close to some I k , i.e. E[$>] = Auk + e 2 , there 
should be a splitting as follows: 

(15) $(t,x) = I k (X(t)x) + {small error} , 

and the goal is to understand the lower dimensional dynamics of the parameter 
X(t). Plugging the ansatz (fT5)) in the action (fTTj) yields the following effective 
Lagrangianj for A(t): 

A 2 

(16) £[*](*) = C k -j{t) + Airk + {small error} , 

A 

where A = 4| and where the normalization constant Ck is given by the explicit 
integral (note that this is only finite for 2 ^ k): 

Ck = - \ I \\rd r I k f dx . 

Here || • || 2 is the norm on M. 3 . The effective dynamics generated by (Q~6|) are now 
given by the formula: 

(17) -*(»-) - {™all error). 

If one were to ignore the contribution on the right hand side of this last equation, 
then the evolution generated by (fll))) would imply that dynamically the soliton ra- 
dius collapses as a linear function of time, or equivalently that A ~ (T** — t)^ 1 for 
some fixed T**. In this sense, the moduli space Mk is said to be incomplete. 

While the above scenario is appealing for its simple geometric and physical mo- 
tivation, it has been rigorously known for some time that it cannot be quite correct. 
This is due to the fundamental regularity results of Shatah and Tahvildar-Zadeh 
[S]-[32] (see also [9]-[10]) which rules out the existence of purely (i.e. linear) self- 
similar collapse: 



3 This is very similar to what is done in the modulational theory of dispersive solitons, and we 
will expound on this in much more detail in the sequel. 

4 It must be kept in mind that this calculation is purely heuristic as, the original Lagrangian 
|(T} itself is only a formal way to derive the equations J2} . 
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Theorem 1.2 (Regularity theory for symmetric wave-maps). Let $ be an equivari- 
ant solution to the equation {2J) with smooth Cauchy data. There exists an e > 
with the following property: Let T* be any time such that this solution is C°° for 
all times ^ t < T* and such that the following condition holds: 

• For any ^ t < T* the energy content J?g/ t ) [$(£)] inside the ball B(t) 
centered at the origin r = 0, is such that lim^T* Ejgu)[<$>(t)] e whenever 
\B(t)\ = (T*-t)-o(l). 

Then the wave-map $ extends past T* as a C°° solution. That is, if the energy 
E[Q(t)\ concentrates at most at a (linear) self-similar rate up to time T* , then the 
solution cannot break down at time T* . 

This theorem shows that the error terms on the right hand side of (|17p cannot 
be ignored, and that any complete theory of how A(t) should evolve must take them 
into account. In fact, the above result leaves the question of breakdown in finite 
time for the equation @ open to a much wider range of possibilities because while 
it gives a necessary lower bound on any possible blowup rate for X(t), it does not 
give any upper bound in case collapse might occur. 

At this point we should further mention that the general wave-map equations 
J3]) have also been studied intensely from an analytic perspective. For the static 
case of ((3|) we point out the references [15] and Chapter 8 of [22] , and the references 
therein. In the case of dynamic solutions and the Cauchy problem, the only general 
understanding of the equations ([3]) that is yet available is for the local theory (see 
[15] ) and the small data scale invariant (global) results of [37] (in the Besov case) 
and [35], [IH], and [33] for the case of small energy. Large data global regularity 
has been conjectured in the case of a hyperbolic space H 2 target, while singularity 
formation has been expected for the 0(3) cr-model for some time. We explain this 
more in a moment. We point out to the reader that this is in stark contrast to 
what is known for the parabolic analog of ([3]) (i.e. the harmonic map heat flow), 
where the global regularity theory at all energy levels is much better understood 
(see [12], 0, 0). 

However, in the case where the dynamic solutions of ([3J possess a large amount 
of symmetry, there has been considerable progress toward our understanding of the 
general Cauchy problem in the case of arbitrarily large initial data. The global 
regularity question was first handled in the work of Shatah/Tahvildar-Zadeh (see 
[31]-[32]) and Christodoulou/Tahvildar-Zadeh (see [9]-[10]), where the context is 
spherical symmetry or more generally k-equi variance. There is also the important 
and closely related work of Struwe (see [33] -[35]), where breakdown is studied in the 
general (symmetric) case including maps into S 2 . From all these works, it is known 
that if the target manifold M. is "geodesically convex" , then symmetric solutions to 
(j3J cannot break down in finite time. Furthermore, this behavior has been shown 
to be stable under small rough perturbations in the recent work of [20] . Finally, it 
is known in general (i.e. without geodesic convexity) that if a symmetric solution 
to ([3]) does break down in finite time, then the singularity formation must be tied 
to the existence of a static solution to ([3]) , and in fact will rescale to a non-trivial 
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harmonic map in the limit. 

While the works mentioned above furnish a great deal of understanding, they 
also leave completely open the issue of whether or not singularities do in fact form 
in the specific case of dynamic solutions to the equation ([2]) . The most convincing 
evidence to date that breakdown does occur in finite time is the analytic work of 
R. Cote on strong asymptotic instability in the energy space, and the many 
numerical studies that have been performed (see for example [2], [T7], [53], [21], 
and [29]). We mention here that the work [23] suggests a universal log- modified 
self-similar behavior similar to ([7]). 

In this work we show that singularities will form in finite time for the critical 
0(3) cr-model in such a way that the stable dynamics is bounded by a log-modified 
self-similar collapse that is not so fai0 from what is predicted by (fTT)) . One of the 
major points of this paper is to uncover the precise analytic mechanism which is 
responsible for this upper bound. 

Before closing this subsection, let us make several remarks. The first is that the 
0(3) cr-model also enjoys many analogies with other more complicated field theo- 
ries such as the (4 + 1) Yang-Mills and the (3 + 1) Yang-Mills-Higgs equations. For 
this reason, the it has been an important testing ground for ideas concerning the 
structural behavior of these more complicated models. We would like mention here 
the work of Bizon, Ovchinnikov, and Sigal for the case of Yang-M ills instantons [3] , 
which proposes a collapse scenario for the critical (4 + l)-dimensional Yang-Mills 
equations similar to what we deal with here. The reader will see that some of our 
methods are inspired by certain calculations performed in that paper. 

Secondly, existence of finite time blow-up solutions had been had been known 
for some time in the case of a super-critical higher dimensional wave map problem 
with Minkowski space R n+1 with n > 2 as a base and a rotationally symmetric 
Riemannian manifold M as a target manifold. The construction of blow-up solu- 
tions is based on existence of fc-equivariant self-similar solutions of finite energy for 
the higher dimensional wave map problem. Such solutions have been exhibited in 
the work of Shatah (see [30]) for the R 3+1 — > S 3 problem. This was later extended 
to other target manifolds in [32] and higher dimensions n > 4 in [BJ. In the latter 
work it was also shown that for n > 7 self-similar blow up can occur even in the 
case when the target manifold is negatively curved. 

Thirdly, an interesting issue that we would like to draw the readers attention 
to here is that the k-equivariant heat flow corresponding to the instance of ([2J we 
study here is known to be globally regular (it is expected that the corresponding 
Schrodinger flow is also globally regular). That is, for the equivariant maps into 
the sphere § 2 with the homotopy index k > 1, the harmonic map heat flow does not 
break down in finite time [14] . The reason why finite time breakdown can occur in 
the wave flow analog of this problem is essentially due to the second order nature 
of the equations. See Remark 1 1 . 71 below for more thorough discussion. 



'In terms of power law type behavior. 
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Lastly, we point out that our work is essentially independent of previous tech- 
niques used for wave-maps, although we find it extremely useful to keep in mind 
that the self-similar blow-up is a priori ruled out by Theorem 1 1.21 However, we do 
refer to that result for the statement of "small energy implies regularity" , which 
underlies much of what we do in the sequel. We again stress that the fundamental 
structure we rely on in this paper is the "quasi-integrable" and "super-symmetric" 
aspects of the static (elliptic) case of the equation j2|). Specifically, the fact that 
such solutions may be constructed by solving the first order Bogomol'nyi equations, 
as opposed to the full second order Euler-Lagrange equations. These aspects enter 
prominently into our analysis of the time- dependent problem. 

In the remainder of this section we will give a detailed discussion of the symme- 
try reduction we use in this work, as well as two separate statements of our main 
theorem. 

Acknowledgments: The authors would like to thank M. Grillakis and J. Shatah for 
valuable discussions of the symmetric wave-map problem, and J. Krieger and W. 
Schlag for valuable discussions regarding their recent work [21j and related stability 
problems of wave equation solitons. We are also deeply indebted to Pierre Raphael 
for pointing out to us a subtle miscalculation of the blowup rate in an earlier version 
of the paper and for other valuable discussions. The authors would like to thank 
the MIT and UCSD mathematics departments for their hospitality while this work 
was being conducted. 

1.1. Symmetric reduction of the problem, and the statement of the main 
theorem. As we have already mentioned, we will restrict our study of the system 
([2]) by enforcing some fairly rigicQ symmetry and "size" assumptions. As seen in 
the introduction, the class of solutions one has access to under these restrictions 
already exhibits some interesting and striking phenomena, and in many ways is 
still quite far from being understood. We now give an alternative derivation of the 
symmetry assumption we use here. This corresponds to solutions behaving rigidly 
with respect to rotations on the base manifold K 2+1 . That is, we require that along 
some fixed time-line (t, 0) € R 2+1 a rotation of 2tt corresponds to a rotation of 2kir 
on the sphere § 2 about some fixed axis. This type of symmetry dictates a more or 
less canonical set of coordinates on the target, which is simply polar coordinates 
centered about the axis of rotation. We write this in the usual way in terms of two 
angles: 

(18) {$\$ 2 } = (0,6»), g = ds 2 = d<p 2 + sm 2 (</>)d<9 2 , 

where we restrict ^ (f> ^ it and ^ 9 < 2ir, with <fi = tt,0 the respective north 
and south poles of the rotation axis. 

With this choice of coordinates, our symmetry assumption boils down to the 
simple relation 9 = fc8 where (t, r, 8) are polar coordinate on the base R 2+1 and 
we mod with respect to 2-7r. In this case, the only remaining degree of freedom is 



"It would be extremely interesting to remove these in some way. For example by either studying 
large initial deviations from an equivariant soliton, or by studying small non-equivariant pertur- 
bations of an equivariant soliton. 
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given by the quantity cj> which can only depend on the variables (t, r). Because all 
of the Christoffel symbols Vfj vanish except in the case i = j = 9, where we have 
by a simple calculation Tg 9 — — \ (sin 2 <fi) , the general system © reduces to the 
single equation: 



(19) 



1 



di<p + {d 2 r + -d r )<j> 



2 $m(24>) 



k e N+ , 



r 2r 2 ' 

where we implicitly enforce the boundary conditions 0(0) = and 4>(oo) = n. 

Before we continue, it is useful for us to record here the formula for the La- 
grangian density ^ under this k-equivariant symmetry reduction and in terms of 
the local coordinates on the sphere (<fi,9): 



(20) 



C[4>] 



-(9^) 2 + (9^) 2 + ^(l-cos(20)) 



In this notation the conserved energy (fT0|) becomes: 



(21) 



m 



(d^) 2 + (^) 2 + ^sin 2 (0) 



rdr . 



The statement of our main theorem is now the following: 



Theorem 1.3 (Finite time energy concentration for wave-maps). Consider the full 
wave-map equation (|6|) with § 2 target under the equivariant restriction to equation 
(|19|) . Then for any integer 4 k, and for any sufficiently small constants < cq <C 1 
with the property that for any e Cq, we can find a set of smooth (in the sense of 
the full map on IR 2+1 ) Cauchy data: 



0(0) 



J ! 



d t 4>(0) = 



o ■ 



with energy size -E[</> e ] = 47rfc + e 2 such that this solution collapses at a finite time 
T** . More specifically, this solution collapses at a rate bounded by a "log-modified 
self-similar" dynamic in the sense that there exists a universal time independent 
profile <j) k , and a real parameter < X(t), such that: 



(22) 



E 



<f>{t 7 r)-0 k (\r) < e 2 
and such that for any < M and times sufficiently close to T** one has the bound: 
M 



(23) 



(T** - t) 



< A(t) < cl 



j V\ ln(T** ~ t)\ 



(T** - t) 



Finally, this type of blowup is stable within the class of initial data in the sense that 
there exists a weighted Sobolev space H s,m (see (|30p for a definition), such that the 
coe ball about (4>q,4>q) in H s ' m also leads to collapse with the same universal profile 
6 and the same bound 
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1.2. The family of static solutions and a modulational version of Theo- 
rem 11.31 As we have mentioned previously, it is well known from work of Struwe 
(see again |33j ) that any blowup of the form described in Theorem 11.31 must in 
fact be a "bubbling off" of a static solution to the equation (fT9|) . That is, after 
rescaling the solution <j>(t) as described in Theorem ll.3[ the resulting profile should 
be a solution to (fl9|) . In the sequel, we will actually take the converse approach and 
give an explicit construction of such bubbling off solutions. This will be done in a 
way which is generally consistent with the decomposition (| 1 5|) of the introduction. 
Our method also naturally shows that this process is reached from a generic (in the 
symmetric sense) set of initial data, and that it enjoys a certain universality which 
is embodied by the blowup rate (f23|) . 



To get things started, we derive the formula for the energy minimizer of the 
(full) action ([2H| . This is just a recalculation of lines (fl"2"|) in the current notation. 
Completing the square in the spatial terms in the energy (|2ip we can write it as: 



E\d>] = TT 



{d t 4>f + {d r ct>--sm{4>)) 2 



r 



(24) 



(a t 0) 2 + (9 r 0--sin(0)) 2 



r 



rdr + 2tt / ksin((j))dr4> dr, 
J 

rdr + Akix 



Thus, one has the universal lower bound Ak-n ^ E[<f\, which can be reached if we 
can find a function I k with the property that I k (0) — and I k (oo) — it, and which 
satisfies the following equations: 

d t I k = , 

(25) rd r I k = fcsin(/ fe ) . 

We shall refer to the solution I k as the harmonic map soliton. A direct calculation 
reveals that the function I k is given by the explicit formula: 

(26) I k {r) = 2tan~ 1 (r' £ ) . 
We also denote: 

(27) I{r) := I k {r) , J(r) := rd r I(r) . 

Note that since the equations (|25[) are homogeneous, the general solution I\ is only 
defined up to a rescaling, as we have already mentioned on line (j 14[) above. Now, in 
terms of these objects we can state the following more technical and precise version 
Theorem II .31 which is what we shall actually prove in the sequel: 



Theorem 1.4 (Modulational version of the main theorem). Consider the reduced 
wave-map equation (|19p with 4 ^ k. Suppose we are given a pair of sufficiently 
small positive constants e, Co with e ^ c\, and an initial data set of the forn^: 

e 

7T 



(28) 0(0) = I + uo, c^(0) = ^ || J\\-J {rdr y J + 9o 



where: 



uq ■ J rdr = , 



^Note that our choice of initial data already requires 2 sC k as for k = 1 we have that 
II (rdr) = 00 • 
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and obeys the smallness condition: 

(29) II K$b) 112pm < c 2 e 2 , 

where we have set: 

(30) 

IIKso)!!^ = £ / {l + r 2 r%dl 90 f + ^ + {d^u,f + ^^]rdr. 

Then we have that the following is true: There exists a continuous time dependent 
parameter X(t) with A(0) = 1, and such that the solution (j) to (|19[) with initial data 
(|28[) splits into the sum: 

(31) 4>{t,r) = I(X(t)r)+u(t,r) , 
where the "remainder" term u satisfies the bounds: 

(32) / [{d t uf + (d r u) 2 + %] rdr < e 2 , 

JM+ r 

for all times the solution exists. Furthermore, there exists a finite time T** such 
that lim t ^T** A(t) — > oo. Finally, this parameter obeys the following bounds for 
times t sufficiently close to T** : 



M I J\\n(T** -t)\ 
(33) t r «S Aft) < c n 4 ^— i 



Remark 1.5. The lower bound in the blowup rate follows easily from the orbital 

stability bound (|32p and Theorem ll.2l Therefore, in the sequel we shall concentrate 

on establishing blowup with the upper bound on line (|33p . The reader should note 

i 

however that the presence of the extra small constant Cq , which may go to zero 
with e, indicates that the true blowup rate is even closer to self similar than the 
y/\ ln(T** — t)\ correction. We will return to this delicate issue in a later work. 



Remark 1.6. The extra decay provided by the bounds (|29p is not essential to what 
we do here and the result also holds in the space H 2 ' . It is assumed here as a 
convenience that will simplify the exposition. However, the extra regularity afforded 
to us in the norm ([30]) will be used in a crucial way. We also remark that the 
norm (|30p is consistent with smoothness of the derivative of the initial data when 
considered as map from R 2 into the pullback bundle U 2;6 R2+i<i>*(T§ 2 ). This is a 
consequence of some simple calculations involving the frame Christoffel symbols 
©• 



Remark 1.7. As we have already mentioned, the blowup mechanism we exhibit here 
is ignited by a spectral phenomenon. The choice of initial data (f!?g|) guarantees that 
the time derivative of the wave-map <po has a "strong" projection onto the "ground 
state" J(r) of the equation (fT9")) linearized around the soliton I(r). This is precisely 
the coefficient in a Riccati equation for the scaling parameter \{t) (see formula 
(f6"3"| below). The Riccati equation generates the first self-similar epoch of collapse 
which lasts on the time interval of size ~ e _1 and, as the projection of the time 
derivative d t (j> on the "ground state" J\ of the modulated soliton I\ goes to zero, 
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is eventually replaced by a more violent accelerated regime leading to the blow-up. 
For this initial phenomenon to take place it is crucial that the linearized ground 
state J is an L 2 (rdr) function and that the projection of the time derivative of <p 
on the ground state J is initially non-trivial. That is, one of the main things which 
makes our analysis possible is that the first ordei@ field quantities (4> — I\, dt4>) can 
not be both orthogonal to the eigenfunction of (|19p linearized around I(Xr) (unless 
one restricts the initial data to a co-dimension one submanifold) . 

In this regard there are some interesting open questions connected with the value 
of the homotopy index. For k = 2 the linearized ground state is still in L 2 , so it 
is likely that an adaptation of our methods is possible. This is important because 
it is this case which is most closely related to Yang-Mills (see the next remark). 
For the unit homotopy class, k = 1, the situation appears to be more complicated. 
In this case the linearized ground state just misses L 2 by a log. The major open 
problem here seems to be whether there is complete instability of the kind stated in 
Theorem 11.31 or if small enough perturbations (in some space) are asymptotically 
stable, with blowup occurring as some kind of "critical phenomena" depending 
delicately on the size of the perturbation. Another interesting thing is that in the 
case of k = 1, there are some numerical simulations which seem to indicate that 
the blowup, while taking place, occurs at an algebraically different rate from (|33[) 
(see again [2]). On the other hand, there are further numerical and heuristic results 
(see [21]) which suggest the validity of the log-modified behavior even in this case 
(k = 1). We believe that both the k = 1, 2 cases of (a possible analog of) Theorem 
1 1.31 deserve further serious investigation in terms of numerics, heuristics, and theory. 



Remark 1.8. Another important issue we call the readers attention to is that in the 
case of the critical Yang-Mills, heuristic arguments as well as numerical evidence 
point to blowup with a modified self-similar asymptotic of the same form as (|33[) (see 
again [3]). In fact, the spherically symmetric reduction of the (4 4- l)-dimensional 
Yang-Mills equations is very closely connected with the homotopy k = 2 case for 
the 0(3) cr-model. This strongly suggests that the methods we develop here will 
transfer to the case of the Yang-Mills model as well, and this will be the subject of 
a forthcoming work of the authors. 



1.3. A few more calculations. Before continuing on, we list here some simple 
formulas involving the unit solitons I k which will be of particular importance to us 
in the sequel: 

O r k i Jlk 

(34) rd r I k = fcsin(^) = k^--^ , cos(/ fe ) = , 

r k „3k -i _ « r 2fe I ik 

(35) sin(2J fc ) = 4^—^ , cos(2/ fc ) = (i - / k)2 . 

Also, in the sequel we will refer to any specific instance of I k as simply /, and we 
remind the reader that we are assuming 4 k. 



^That is, in this notation the second order wave equation (1191 can be written as a first order 
system. 
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1.4. Vanishing of the wave-map. We end this section by recording and proving 
a simple geometric lemma which will be of central importance to us throughout the 
sequel. We will show that a wave-map $ together with its derivative vanish at the 
origin r — when computed in the pair of local coordinates (fT8|) and (r, 9). 

Lemma 1.9 (Admissibility condition for the wave-map <£). Let $ be a smooth 
k-equivariant function from the plane M. 2 into the sphere § 2 . Then if2^k one has 
that: 

(36) \d T <t>\ ^ C^r , < r < 1 . 



Remark 1.10. Note that the condition 2 ^ k for estimate (|56"|) is crucial, as the 
formula (|26[) shows for unit homotopy class instanton I 1 . 



Proof of the estimate (J36J). Our first step is to establish that d r <f> is continuous and 
vanishes at r = 0. 

First of all, notice that along any fixed radial line = const, the vector-field d r 
is a continuous section of TM. 2 . The same is true of the field ^Oq. Furthermore, 
one has that: 

lim d r = — lira — d<$ . 

e=o e=f r 



r— >0 



Therefore, by continuity we must have that: 



lim \\d<5>(d r )\\ 2 = lim \\d^(-d e )\\ 2 



+0 

12 



where || • || 2 refers to the metric (fT5j) . Computing both sides of this last equation, 
we see that not only is <9 r </> continuous (and hence bounded) on the interval [0, 1], 
but that we also have: 

in ,, -,. fc|sin(</>)| 
(37) lim \d r <t>\ = lim -! . 

Using now the fact that (f>(0) — to write <fi(r) = f£ d y cf>(y) dy, upon substitution of 
this integral into the right hand side of (j37| , we see from the fundamental theorem 
of calculus and the condition 2 ^ k, that we must in fact have d r (f>(Q) = 0. 

It remains to show that d r <p vanishes uniformly (with non-uniform constant) in 
r. To do this, we compute the Dirichlet energy: 



1 



e(0) = gij d r &d r & + - s d e &d 9 & = \d r <p\ 2 + — |sin(0)| 



k 



2 



2 



r 



This is a C°° function on R 2 which depends on the radial variable only. Further- 
more, we have that e(</>)(0) = 0. Thus it is necessary that: 

\e(4>)\ < Cr 2 , < r < 1 . 

for some constant that depends on <j>. In particular, we have the bound ([5t?)) . □ 
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2. Some notational conventions and an overview 

In this section, we will first list some standard notational conventions that will 
be useful throughout the sequel. We then give a quick technical overview of the 
main result. 



2.1. Some notation. Throughout this paper, we shall employ the standard no- 
tation A < B to mean A ^ CB for two quantities A and B, where C is a fixed 
constant. There is no uniformity in this notation for separate instances. That is, 
separate occurrences of < on the same page will not necessarily imply C is the same 
for each. Another, less standard, notation which will be of great use is the following: 



Notation 2.1. For any pair of non-negative integers 0^m,nwe will denote by 
pm,n an y £7°° function which satisfies the following bounds: 

Urd r yF m ' n \ < C, 



(1 + r) m +™ 

We will also use a shorthand notation for the case m = 0. Here we shall set 
F n = F°< n , so that we have: 



\{rd r ) l F r ' 



< 



(1 + r) n 



We also denote the A rescaling of these functions by F x ' (r) = F m ' n (Xr), and 
similarly for F£. Finally, let us remark that different instances of F ™ : ™ on any line, 
or between lines, can mean separate functions. 



This notation will occur so frequently in the sequel that is will be useful for us 
to record here several instances which involve either time or space derivatives, or 
multiplication by powers of r. Collectively these are the following, where we assume 
j G Z is such that —m ^ j in the first identity and 1 m in the third: 

(38) 

All of these are immediate from Definition 12.11 above. 

Also, in the sequel we will in general use the A-subscript notation to denote the 
A rescaling of a given function. For example I\(r) = -f(Ar). 



2.2. An overview. As is perhaps already clear at this point, our method for es- 
tablishing Theorem 11.41 is to control directly a certain modulational equation for 
the time dependent scaling parameter A(t), and to show that this evolves according 
to a blow-up ODE. Thus, in this sense our work is closely related in spirit to the 
modulational stability approach originally pioneered by M. Weinstein (see [31]) and 
later sharpened by Buslaev and Perelman (see [S]) to study solitons dynamics of the 
focussing non-linear Schrodinger equation. The major difference however is that we 
are actually trying to show that there is an extremely strong asymptotic instability. 
This of course introduces a serious problem when one tries to control the non- linear 
equation (|19| linearized around the modulated soliton I\ . Nonetheless we begin by 
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using the decomposition <)) = I\+u from line ([31]) , and then linearizing (|i9|) around 
h: 

(39) d 2 u + H x u = -I x + JV(u) , 

where the Hamiltonian is given by H\ = —d 2 — r~ 1 d r + Q\{r), and the nonlinear 
term Af(u) containing quadratic and higher order terms in u (also containing a 
factor of r~ 2 ). 

Our first task, dealt with in Section |4] is to prove orbital stability of the modu- 
lated soliton I\ under the condition that the "radiation" part of the solution (i.e. u) 
is orthogonal to the function J\ = rd r I\, which is the unique eigenfunction of the 
Hamiltonian H\ . The latter is a consequence of the fact that I\ realizes an absolute 
minimum of the energy (|2ip associated with the full nonlinear problem (| 1 9[) . The 
orthogonality condition provides us with an ODE for the scaling parameter X(t), 
which is coupled to the radiation term u: 

(40) x(^(J,J) - (u(X~ 1 r),rd r J)^ = (A(0 t ,J A ))A 2 , 

while the orbital stability statement will give us a very weak control of the remainder 



u: 

r r ;.2 i 

2 



(41) 



(Mf + (d r uf + ^U 2 



rdr < e 



Due to the coupling between the scaling parameter X(t) and the radiation u, to 
control A(t) to the extent that we can show A(i) — > oo in finite time requires much 
better control on the radiation term u. The usual procedure for dealing with this is 
to scale out the modulational parameter \{t) at each fixed time, so time independent 
spectral methods can be used to control the linearized equation. This procedure 
works well if one can prove that there is a slow limit of the parameter X(t), but it 
obviously causes a catastrophe if X(t) grows rapidly. In this truly non-linear 

approach is needed. 

The reason why standard non-linear estimates, for example the kind used to 
prove orbital (Lyapunov) stability (i.e. (|41 [) above ), are not sufficient to reach the 
blow-up time ~ e _1 is basically due to the fact that their application to the ODE 
(|40p is not truly scale invariant. That is, the use of fixed time estimates which result 
from orbital stability analysis causes a loss (of scaling) when one integrates over 
time. Such integrations seem unavoidable when analyzing (|40[) . To overcome this 
problem requires uncovering a non-linear dispersion phenomenon^ in the equation 
(|39| for the radiation term u. That such a dispersive process indeed takes place is 
in some sense the miracle of the equation (|39j) . More specifically, as the soliton 7a 
collapses it actually repels the excess radiation away from the origin. This evacu- 
ation process only causes the soliton to collapse at a faster rate, and it is what is 
ultimately responsible for the acceleration of self-similar behavior governed by the 
LHS of the blow-up rate {33]). 



9 A non-linear dispersion phenomenon (of a different nature) has been observed and used by 
Merle and Raphael in their work 28 on the blow-up analysis for the critical focussing non-linear 
Schrodinger equation. 
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One of the most interesting issues in this paper is the mechanism by which this 
"repulsive" behavior of the linearized equations manifests itself mathematically. 
This is where the "quasi-integrable" system aspect of the static version of (fT5|) 
comes in. As we have already discussed in the introduction, static solutions to (|19[) 
are generated by the first order Bogomol'nyi equation (|25p . When one linearizes 
(fT9|) around these static solutions, the corresponding Hamiltonian H x splits as a 
product of two first order operators which are adjoints of each other. That is: 

Hx = A* X A X , 

where A x is the linearization of (f25|) . The reason why this splitting is so useful is 
that the Hamiltonian H x also possesses its super-symmetric companion: 

Hx = A X A* X , 

with A* x and A\ being the analogs of the creation and annihilation operators, and 
H\ and H\ related to each other according to the remarkable intertwining relation: 

(42) A X H X = H X A X . 

Such a splitting elucidates the non- negativity of H x , and identifies the function J x , 
which is the kernel of A Xl as the ground (vacuum) state of H x . In addition, the 
intertwining property (|42[) allows us to simply conjugate the problem (|39[) to one 
whose linear part involves the more manageable Hamiltonian Hy. 

(43) df(A x u) + H x (A x u) = -A x (I x )+A x Af(u) + [df,A x ]u. 

The Hamiltonian H x , which is obtained from H x by the process of "removing" its 
ground state, is of the explicit form H\ = —d 2 — r~ 1 d r + V\(r) and involves a 
space-time repulsive time-dependent potential. This means that for the problem 
(|43|) one may proceed via purely physical space methods, and there is no difficulty 
in handling extremely violent growth of the scaling parameter \{t). What we can 
do is to establish quite strong (i.e. scale invariant) integrated and fixed time energy 
estimates (i.e. so called Morawetz type estimates), while keeping precise track of 
the influence of the source terms on the right hand side of (|43| involving the scaling 
parameter A. To undo the conjugation procedure embodied in (|43|) . we only need 
to use the fact that u is orthogonal to the kernel of A x , because through a little 
elementary functional analysis this allows one to turn our Morawetz estimates into 
ones involving only the term u (as opposed to A x u). Once these estimates are 
established it is possible to show, through a somewhat lengthy calculation, that 
after a long self-similar epoch where CoA ~ eoA 2 , the modulation ODE for \(t) 
enters another monotonic regime where it takes the final form: 

C7 A(i) ~ e X 2 (t)-\ 2 (t) J o(^(8))ds . 

It is this ODE which leads to the blow-up, and also gives the tight upper bounds 
in . The reader should compare this last formula to the blowup ODE for mod- 
ulated Yang-Mills instantons derived in [3] through heuristic arguments. 

We now turn to the details of all of this. As is common in this type of work, many 
of our assumptions will be bootstrapped. We shall follow the outline: 
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• In Section [3] we discuss the Hamiltonians H\ and H\. 

• In Section 2] we derive the basic ODE for the scaling parameter \(t) and 
prove orbital stability statement. 

• In Section [5l assuming certain estimates on the "radiation" part of the 
solution, we obtain a much refined closed form (i.e. without explicit de- 
pendence on the u) of the modulation ODE, establish its monotonic and 
algebraic properties, and prove blow-up along with an explicit rate bound. 

• In Section [6j assuming the monotonic properties of the scaling parameter 
X(t), we prove the integrated space-time and fixed time bounds for the ra- 
diation by making use of the conjugated Hamiltonian H\. 

• In Appendix lAl give some further explicit computations needed in the anal- 
ysis of the blowup ODE. 

• In Appendix [B] we establish some general coercive properties for the class 
of first order operators related to A\. 



3. The Linearized Equations and a Basic Spectral Calculation 

Our purpose here is to derive and record certain calculations involving the lin- 
earization of the equation (|19|) around a time dependent modulation of the soliton 
I(\(t)r). That is, we decompose the full solution as: 

(44) <j>(t,r) = I(X(t)r)+u(t,r) . 

This yields the following set of formulas for equation (|19p linearized around I\ : 

(45) d\u + H x u = -h + N(u) , 

where we have set: 

,, , fc 2 sin(2J A ) , , NX fc 2 cos(2/ A ) , 1 , NN 

Here H\ is the linearized Hamiltonian: 

H x = -a 2 -i9 r + ^cos(2/ A ) , 

(46) = A* X A\ , 

where the first order operators A, A* are given by: 

(47) A x = -d r + - cos(J A ) , A* x - d r + - + - cos(7 A ) . 

r r r 

The spectrum of Hx, defined as a self-adjoint operator on L 2 (rdr) and obtained 
by taking the closure of C^°(R + ) in the graph norm of H x , is easily computed via 
the factorization (|46|) , or via the knowledge that the ground state Ix of the (static 
form of the) equation (|19p is unique. There is a unique eigenfunction, which has 
zero energy, and it is given by the formula: 

(48) Jx = rdrlx = fcsin(7 A ) . 
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In particular, J x solves the first order "linearized BogomoPnyi equation" : 

(49) A X J X = 0. 

The absolutely continuous spectrum fills the half-line [0, oo). 

Of primary importance for use here will be the conjugate operator formed by 
A\A* Xl the super-symmetric companion H x of H x . 

~ 1 k 2 + 1 2k 

(50) H x = A X A* X = -d 2 - -d r + — ±— + cos(/a) , 

= -d 2 r -^d r + V x (r) . 

Recall that the Hamiltonians H x and H x are related via an intertwining relation 
(l42f . As opposed to H x , the spectrum of H x has only an absolutely continuous 
component filling [0, oo). What is more important is that while in our application 
the Hamiltonian H x will be time dependent, it has a remarkable structure which 
allows one to prove strong local energy decay estimates even if the parameter X(t) 
grows in an unconstrained fashion. The needed properties follow from entirely 
elementary calculations and are as follows: 

(51) V x ^ (fc ~ 1)2 , (Positive) , 

4fc 2k 2 
d r V x = — ' ^ + cos(/ A ) H s- sin 2 (/ A ) , 







2(/e 2 H 


-1) 






2(k- 


I) 2 





(52) ^ , (Space-Repulsive) 

A 2k 2 

(53) — d t V x = — ■ —5- sin 2 (/a) , (Time- Repulsive) 

A r l 



We conclude this section by refining the decomposition (|44|) . This will be ex- 
tremely important for us in the sequel, and it will ultimately lead us to the precise 
asymptotic (|33[) . What we will need to do is further decompose the radiation term 
as: 

u = wo + w , 

where the "leading term" w is meant solely to eliminate the term I x on the right 
hand side of (|4"S"|) , while at the same retaining the orthogonality relation: 

(54) (w ,Jx) = 0. 

The equation we use to generate wq is (the reason for this will become apparent in 
the sequel): 

(55) A x (H x w ) =H x (A x w ) = -A X {I X ) . 

To further the computation, we use the fact that A X J X — to write: 

Aa(Ja) - jA x d t {rd r h) = ^A x rd r (J x ) = ^A X (J + rd r J) x . 



Therefore, by peeling off the common factor of A x from both sides of (|55|) it suffices 
to solve: 

A 2 

H x w Q = - — (J + rd r J) x . 
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We now use the ansatz wq = A 2 A~ 4 [-f£T(Ar) — 7 Ja] , where 7 is a normalization 
constant which will be chosen as to satisfy (|54p . Eliminating common factors, and 
rescaling the spatial variable we find that: 

Hx(K) = -(J + rdrJ) , 

from which a direct computation shows that K(r) — ^J{r) is the desired solution. 
Therefore we have that: 

(56) w {t,r) = ^(t)-(aJ x (r) + b(r 2 J) x ) , 
where the coefficients are given by the explicit formulas: 

(57) a = -i(J, (r V)).||J|| Z2 2 Mr) , b=\. 



Before ending this section, let us translate the abstract function notation on line 
(|38|) into some specific bounds which will be used many times in the sequel. These 
are a consequence of simple explicit formulas, and the restriction 4 ^ k: 

(58a) w = A 2 A- 4 (F A 4 ' 4 + if 2 ) , d t (w ) = X\\-\F A X A + F^ 2 ) 

+ A 3 A- 5 (F A 4 ' 4 + F A 6 ' 2 ), 

(58b) |J A | < F 4 X 4 , \8 t J x \ < AA- 1 ^ 4 , 

\dUA < (lAIA-^A- 2 )^ 4 , 

(58c) 

d t (Ax) = A^ 7 ' 9 , \d 2 (A x )\ < (lAI + A^- 1 )^ 9 , 

MK' n ) = AF™ _1,,1+1 . 
We note here that we are assuming 1 ^ m for the last identity on line (|58c|) . 



4. Orbital Stability 

We now begin with the first step in our proof of Theorem 11.41 This is to show 
that one can make a rough decomposition of the full wave-map <p into a bulk piece 
which is a rescaled soliton, plus a small remainder which we can estimate in a cer- 
tain energy space. Of particular importance to us will be that we can construct 
this decomposition in such a way as to retain a certain orthogonality between the 
bulk piece and the small "radiation term" . Of course, this is precisely the modu- 
lational approach to orbital stability first pioneered by M. Weinstein in his study 
of the non- linear Schrodinger equation (see [33] )• What we intend to prove is the 
following: 



Lemma 4.1 (Orbital stability with orthogonal decomposition). Suppose that (j) is 
a solution to the problem (|19p. and suppose that initially the Cauchy data for (f> 
decomposes as: 



(59) 



<M°) = ^Ao + Uo , 



cW(0) = 0o , 
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Assume that the energy satisfies E[(f>] = Ank + e 2 , with e chosen small enough. Then 
as long as the solution 4> exists there is a time dependent parameter < X(t) < oo, 
with A(0) = Ao, and the property that the following conditions hold for all times of 
existence provided that they hold initially: 



EM = I 



{dt4>f 



(d r uf + ^U 2 



rdr < e 



(60) 

(61) = (u(t),J m ) . 

Here u is defined by the relation u = <j> — I\. In addition we have that: 

l A 2 



(62) 



< 



Finally, we remark that (conversely) the full wave map is C°° up to any time T 
as long as X(t) < oo and (|60|) holds true for any t £ [0, T]. 



Proof of Lemma The proof essentially reduces to denning an appropriate equa- 
tion for the evolution of X(t), basic existence and uniqueness theory of ODEs, fol- 
lowed by the coercive estimate (|162|) proved in Appendix [5] We note here that the 
last remark of Lemma 14.11 follows from the local "small energy implies regularity" 
statement for symmetric wave- maps, and is for instance contained in Theorem 1 1.21 

We will use here the following observation: Notice that if A — > oo or A — > 
on some time interval such that (pO)) holds, then the full wave-map $ must break 
down also on that time interval because its energy concentrates at r = or r = oo 
respectively (in fact, by finite speed of propagation it cannot happen that A — > 
in finite time, so any "blow-up" of this type must occur at t = oo). 

Now, for a strictly positive real valued function of time A(i), we define the 
equation: 

(63) A(2(/,J) + (0(A- 1 r),r9 r J)) =-(0 f ,J A )A 3 . 

Notice that as long as the wave-map <f> exists and is smooth, this equation is of 
the form a(X,t)X = (3(X,t) for two C 1 functions (X,t) of the variables a, [3. We 
now construct X(t) from (|63[) via a simple bootstrapping procedure. Our goal is to 
provide a strict lower bound for |a| so that we may simply apply the usual existence 
theory to (|63| which then produces X(t). 

First of all, notice that a simple calculation involving the identity: 

(I, rd r J) = -2(1, J) - (J, J) , 

shows that the bound ((60|) initially implies that: 

\a(X ,0) + (J,J)\ < Ce, 

for some, possibly large, constant C which we choose in a moment. Therefore, by 
continuity there exists a small time T* such that the solution <j> exists and is regular 
on [0, T*], and one has that a solution X(t) to equation (|63|) exists, is contained in 
(0, oo), and obeys the bounds: 

(64) \a(X(t),t) + (J, J) I 2C*e . 
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The heart of the matter is now the following: We will show that if T* is any time 
such that the above holds (i.e. existence for (j> and A and the bound |64|) ). then 
we must also necessarily have the conditions (|6ip and (|60[) , as well as the following 
improvement of (|64[) : 



(65) \a(X(t),t) + (J, J) | Ce . 

By continuing this process, we will have shown that on any time interval [0, T*\ 
such that the solution <fi exists and is regular, there is a continuous solution of (|63|) 
X(t) € (0, oo) such that A(0) = Ao and (f6"Tj) - ([6T))) holds (assuming, of course, that 
these conditions hold initially). 



We now show that existence and (|64[) , implies ([61]) and ([60]) , and that these two 
together imply the improved bound (|65p . Everything rests on the orthogonality 
condition (f6Tj) . We define u — <fi — I\ and use a few simple integration by parts to 
compute that: 

4:( u i J \) = ( u t,J\) + ^;(u,(rd r J)x) 1 

at A 

= A( 2 </,J) + (^A-V),r^J>) + (0 t ,J A ) , 
= , 

where the last line follows from the assumption that A solves (|6"5)l . Therefore, since 
(j6Tj) holds initially, it holds on [0,T*]. 

It remains to show (I5D]) . and that this bound implies (f6"S"]) . In fact, given the 
form of a this latter implication is immediate from (|60p and the Cauchy-Schwartz 
inequality, where the sufficiently large constant C is chosen according to the im- 
plicit constant appearing on line (|60[) . 

Thus, we have reduced things to proving that, assuming that orthogonality con- 
dition (|6ip holds, (|60p a lso holds with a fixed implicit constant which does not 
depend on the size of A or C. This is where the condition E[(j>] — E[I\] = e 2 enters. 
Notice that this equality follows from the conservation of energy and our assump- 
tion E[<fi] = Auk + e 2 . Computing the energy difference and using line we have 
that: 



c 2 



{d t 4>) 2 + (d r (j) — — sin(0)) 2 rdr 
(d t 4>) 2 + (d r u - - cos(Ix)u) 2 



rdr — tt1Z\(u) 



where we are defining higher order integrated non-linearity: 

f k 2 f k 2 

7Z\(u) = / {d r <p sin(0)) rdr — I (d r u cos(I\)u) rdr . 

Jr+ r ' Jr+ r 

Now, using the first coercive bound (|162|) of Appendix [Bl we have the estimate: 
(66) E Q [u] < e 2 + |^ A ( U )| . 

Finally, using the simple algebraic formula for the difference of squares, the equation 
(|25]l . and writing the first few terms in the Taylor series for sin(/A +u), we easily 
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have the nonlinear bound: 

I^a(«)| < / i\d r u\ + M) • (-) rdr , 

(67) < (S M) S , 

where the last line follows from Cauchy- Schwartz and the Poincare type estimate: 

(68) \u(r)\ 2 < 2 ^ (<9 r u) 2 rdr^) * • f| rdr 

The bounds (|66 p - (p71) taken together show that we may conclude (f5U|) for e small 
enough and some universal implicit constant (e.g. this can be shown through an- 
other continuity argument). 

Having established ([60]) and (|65|) we can easily show from ([63]) that: 

l-l < < 
' A 2 ' ~ 

This concludes our demonstration of Lemma 14.11 □ 




5. The effective evolution and the main blowup argument 

We now begin in earnest the proof of the main Theorem (|1.4p . This centers 
around computing a more effective form of the modulation ODE (|63|) . This will be 
followed by an ODE analysis giving the desired blow-up together with its asymp- 
totic profile. The main technical result of this section is the following: 

Proposition 5.1 (Refined structure for the modulation equation ([63])). Consider 
the scaling parameter X(t) which is defined through Lemma \4-1\ and equation (|63|. 
Suppose that the initial data (|59p is given according to Theorem \l-4\ Then on a 
time interval where X(t) £ (0, oo) and t € [0, e~ 4 ] the scaling parameter A(i) satisfies 
a first order ODE: 

(69) (C -ei(t))A(t) = e A 2 (i) - X 2 (t) f S(s) ds . 

Jo 

Here Co = (J, J), and cq, ti{t) and £(t) obey the conditions (eo is fixed): 

(70) |eo-^| < c e , 

(71) |ei(?)| < e , 

(72) \£(t)\ < cqc 2 + 4 sup A 4 A~ 7 (s) . 

Here e, cq are the small constants from line (|29p. In addition, one has the following 
"structure bounds" for the acceleration of X(t): 

(73) |A(f)-2^(i)| < Ccl^(t) + C(c e 2 + sup ^(s))X 2 . 

A A O^s^t A' 

We remark that the constant C is universal, i.e. independent of cq and e. 
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Remark 5.2. The reason the explicit constant C appears in the estimate (|73|) is 
merely a notational convenience. The bound (|73|) will be one of our main boot- 
strapping assumptions in the sequel. 

The remainder of this section is divided into two parts. First, we will show that 
the identities and estimates (jSilf - flTB")) imply that the parameter X(t) goes to infinity 
at some time T** £ [0, e~ 4 ] as long as e is chosen small enough. We then establish 
the bounds (|33|) on the asymptotic rate of X(t) as t — > T**. Finally, in the last 
subsection we state the "Main Estimate" (a certain fixed time energy estimate) of 
our paper, and we use it to derive all of the assumptions ([55]) - ([73")) . This main 
technical estimate will be the subject of the final section of the paper. 



5.1. Proof of the blowup and the universal bound for A(i). Using the equa- 
tion (|69|) . as well as the assumptions (f70 |) ~([73 |) we now show that A must blowup 
in finite time. The basic idea is the following: without the contribution of the in- 
tegral on the right hand side of (]69[) the desired blow-up would occur in finite time 
~ Co^q 1 in a self-similar Riccati fashion. Therefore, the only problem is that one 
must guarantee the integral term (which adds a negativd 1 "*! contribution, creating 
a delay effect) does not interfere to the extent that A is driven to zero too quickly 
before the blowup can occur. The fact that a priori (i.e. again by Theorem II. 2\i 
at any supposed blow-up time one must have that AA~ 2 — > 0, indicates a very 
delicate balancing between the two main terms on right hand side of ([M]) . This 
constitutes one of the main technical difficulties in this paper, and why many of 
the estimates which appear in the sequel are so involved. To establish the needed 
control, we will first show that the initial self-similar behavior, approximated by the 
ODE Co A = eoA 2 , forces X(t) into a different monotonic regime, where in particu- 
lar the terms A 4 A -7 dominates the error estimate for £ (i.e. ([72]) ). At that point 
blow-up is assured. Our last task is then to analyze the balance of the integral and 
eo terms on the right hand side of (|69[) . and to derive the precise blow-up bounds 
(|3"3")l from this. Along the way, we will show that all of this can be accomplished 
before the time interval [0, e~ 4 ] expires, so that we still have access to all of the 
structure included Proposition 1 5. II 

We now proceed with the details outlined above. The first main thing is control 
the size of the interval where one cannot guarantee good bounds on £ (s). The key 
to this is to show that A 4 A~ 7 becomes monotonic soon enough and with enough 
force to cover the constant cqe 2 which is lost on line (|T3"|) . Luckily it is not hard to 
show that these two things happen at essentially the same time. Computing the 
time derivative and then using the identity (|73[) we have: 

d ( A 4 \ A 5 ^ 2 A 4 \ A 3 

(74) Jt{y) = A^ + °( C °V + Co£ + V (S) ) A^ ' 

We now let To be the first time such that: 

(75) A 4 A- 7 (T ) = Cc\e 2 , 



That the contribution is overall a negative one follows from Theorem 11.21 
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for some large constant C which is larger than twice the implicit constants in (|72|) . 
First we argue that such a time must occur. 

If such a time does not occur, then if cq is chosen small enough (we remind the 
reader that this is done by simply choosing initial data according to ([29]) , and does 
not affect the size of constants in estimates like (|72p ) one sees immediately from 
([72| and (|69|) that the following bound holds for all times: 

(76) ^r(i) > ^ £ o- / Ccoelds , 



o 



where C is fixed and independent of co,e. A simple argument, which we leave to 
the reader, shows that if A(0) = 1 then this last inequality implies that A — * oo at 
some finite time T** ^ 4Coe n ~ 1 and that in addition: 




— it) > 

for all t G [0, T**). In that case however we would also have that: 

A( t )^ X( t ) > (j-) 4 AW - oo as t^T**, 

which shows that in fact the time T defined above must occur. What's more, a 
simple analysis of the previous argument shows that this time must also satisfy the 
conditions: 

(77) T < 4C e^ , AA- 2 (T ) ^ -^e . 

4C<o 

Notice it is also clear that for all times t € [0, T ] we have X(t) > 0. 

Now, by a direct application of the orbital stability bound (|62p we have that for 
any time where the solution exists there is the bound: 

Applying this on the right hand side of f74|) we see that for all t To, and as long 
as A 4 A _7 (£) ^ Ccqe 2 , we have that: 



By bootstrapping this argument, we see that A 4 A 7 is monotonically increasing for 
all times t T . 

Before continuing, with the proof of blowup, we pause for a moment to upgrade 
the bound ([75]) . This will be used in a crucial way in the sequel (see the proof of 
Proposition 16.31 in Section [6]). We claim that there exists a time T* = To + 0(1) 
such that the following improvement of (|75|) is valid: 

(78) A 4 A~ 7 (T*) = c e . 
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Again by contradiction, if this were not the case by equation (|69|) . the definition 
(f75l) of To, and the time bounds (f77l) . we would have a bound of the form: 



A , , 1 /"* 1 



for times T* = To + 0(1) as long as Co <C 1. Integrating this, and again applying 
(|75[) as well as the orbital stability bounds (|62)l to the term A _1 (To), we arrive at 
the inequality: 

Using the condition that e ^ Cj we see that such a bound must expire in t — Tq = 
O(l) time. Finally, notice that by using the time bounds {77]), the definition of 
T ([7SD, the definition of T* {75]), as well as the relation T* = T + 0(1), we may 
integrate the quantity A 4 A 7 over [0, T*] to obtain: 

r T * A 4 

(79) / ^(s)ds sCe. 



We now return to the main thread of the blowup argument. So far we have 
achieved the following. There exists a time To such that for To ^ t we have: 

• A 4 A~ 7 (t) is a monotonically increasing function. 

• A 4 A~ 7 (<) ^ CcqE 2 , where C is at least twice the implicit constant in the 
bound (1721). 



• The time To is associated with the bounds ([77|) . 

We are now at the point where blowup with the rate bounds (|33[) is assured. From 
the above conditions, we have that for all times t ^ To: 

(80) £ < c|a 4 A~ 7 . 

Writing: 

j(t) = -co + / £(«) ds , 
Jo 

we have that AA~ 2 ~ —7. Differentiating 7 and using the bound ([50)1 we see that: 

7 < Po> - -7%V 

By dividing through by — 7 3 and integrating both sides of this last inequality over 
the interval [To, t{\ we arrive at the estimate (we may assume that —7 through- 
out this argument, as will become apparent on the next line): 

' < -J— + cJIn(A(ti))-4ln(A(T )) , 



7 2 (ii) ~ 7 2 (r ) 

< e- 2 + C |ln(A(i 1 )) 
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where to obtain the last inequality we've used the second bound on line 1TT[) as well 
as the fact that 1 ^ A (To). We now recast this last expression in terms of A: 

(8i) i | < ■ 

V e - 2 + c "ln(A(ti)) 

Integrating this last line over time intervals past [0, To], we see that within 0(e _1 ) 
time we must have e -2 ^ Cq In (A(t)). Therefore we are assured of the bound: 

x 2 - s (t) < A ' (t) < X(t) , < S < 1 , 



inside of some interval [0, Ce" 1 ], for a uniform constant C. This is a Ricatti type 
inequality, which easily implies that X(t) — * oo in 0(1) time starting with 1 A at 
the first time where it holds. Let us call the blowup time T**. 



Finally, we need to recover the rate bounds ([33]) . By Remark 11.51 we need only 
establish the upper bound. From ((81]) we have the inequality: 

for times t sufficiently close to T**. Making the substitution A(i) = e s and inte- 
grating from t to the blowup time we have the estimate: 

i r°° 
Co*(T**-t) < 2 / s 2 e~ s2 ds , 



ds 



VMAW) , o(i) 



< 



A(t) A(<) v /ln(A(t)) ' 



\(t) ■ 

where the last line above follows from the well known asymptotics of the error 
function. The above identity easily implies that as t — *■ T** one has the bound: 



A(i) < d ^W^W\ 



This establishes the upper bound in (|33j) . 



5.2. Derivation of the Main ODE and its Structure. We now derive the ODE 
(155]) . as well as all of the accompanying structural assumptions ([TO]) - ([73"]) . This will 
require a certain fixed time energy estimate which is the main technical estimate 
of the paper and will be proved in the following section. We start by recomputing 
the modulation equation ([55]) using the splitting <fi = I\ + u. Differentiating the 
orthogonality relation ([t?T]) with respect to time, we have the simple identity: 

(82) (j A ,J A ) = (4>t,J x ) + (u,jx) ■ 
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Differentiating one more time, and rearranging things with a little help from the 
equation cf>u + H\u — M(u), we have that: 

(83) (h, J A ) = 2(d t u, j A ) + (u, J A ) + (Af(u), J A ) , 

(84) = 2d t {u, j A ) - (it, J x ) + (Af(u), J A ) , 

where Af(u) denotes the nonlinearity: 



/c 2 sin(2/ A ) , , . , fc 2 cos(2J A ) j 
N(u) = . (l _ cos(2u)) + J-^Z.( u --sin(2 U )) 



(85) = ^— ^ -u 2 +N{u) . 

Notice that: 

(86) \N{u)\ < . 

Next, a short computation shows that we have the identity (recall the definition of 
J A from line ([27)) ): 

<'»•*> " (i-4)< j - j ^ c 4(f) a " ■ 

where we have set (J, J) = Co- Therefore, from these last two lines as well as the 
identities (|83|) — (|84[) . we have our two main structural equations: 

A 2 / \ 

(87) A — 2— = [2(d t u, j A ) + (u, Jx) + (Af(u), J A )J A 3 , 

(88) C Q j t (J^J = 2d t [(u, j A )A] - 2(u, j A )A - (u, J A )A + (M(u), J A )A . 

The first equation ([87)1 is sufficient for us to prove the bound ([73)) . The second 
equation ([88)) will yield (I69[) upon integration. Doing this over a time interval [0, t] 
we have the identity (recall that A(0) = 1): 

(89) C ±(t) - 2{u(t)J Ht) )\(t) 

ds , 



X(0)(C -2(u(0),rdrJ)) - J (<- _( s ) + £( s )) 



where we define the constant C* as follows: 

(90) C, = - fc 2 ( ^ + y j) ^ 2 ,sin(2/ A) .J A 

+ A 2 ( (aJx + b(r 2 J)x) , (rd r J)x) - \ 2 (rd r (aJ A + b(r 2 J)x) , (rc> r J) A ) . 

In Appendix A, it will be shown that C» = 0. This indicates that the precise rate 
of blowup in the inequality (|33| is quite delicate. We'll return to this in a later work. 
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The error term £(s) is given by the expression: 

(91) £ = 2{w,J x )X + (w,J x )X + (w ,(\-2^-)(rd r J) x ) 

( W -(2u,o + w ) ;Sin(2/A) , Ja)a _ ( _v" (u);Ja)a . 

Here the terms w ,ui refer to the decomposition on line (|56p above. We list this 
here again for the convenience of the reader: 

A 2 

(92) u = -^(aJ x + b(r 2 J) x )+w= w + w , 

where the constants a, b are derived on lines (|57[) . Notice that the constant C» 
defined by (|90[) arises from the expression (and a few integrations by parts): 

C*\ 4 \- 7 = 2(w ,j x )\ + (w ,J x )\ 

- (^i,sin(27 A )- J X )X - (w ,(X-2\)(rd r J) x ) . 



Before commencing with the proof of the estimates (|72j) and ([73)1 , let us first 
derive from (|89[) the identity (|69p. and also the conditions (|70p - (l7T|) . First of all, 
notice that from the orbital stability bounds (p0|) we have that: 

\(u,J X )X\ < j \\r^ 1 u\\ L 2 {rdr) ■ \\ (r 2 d r J) x \\mrdr) , 

< e A 

~ A 2 

Therefore, we may define ei on the left hand side of line (frJ9")l as ei = 2A 3 A~ 1 (w, J x ) 
and we have (I7T1). 



Similarly, at the initial time, the identity ()82jl gives the relation (we are assuming 
A(0) = 1): 

\{0)(J,J) = ((f>t(0), J) + X(0){u(0),rd r J) . 
Substituting into this last relation the initial data (|28p . and using the smallness 
condition (f2"9"| we easily have that: 

\(0)(C -2(u{0),rdrJ)) = - + O(c e) , 

which gives the condition (|70|) . Finally, notice that this last line also implies the 
initial expansion: 

A(0) = -j^jr + Oicoe) . 
ir (J, J) 

Plugging this into the first term on the RHS of formula (I57| , and using the bounds 
(f2"5)) on our chosen initial data (|2"8|) to estimate the remaining terms, we see that 
(fT3")l holds for the initial time t = 0. 
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It remains for us to derive the bounds (|72|) and (|73|) from the identities (|9i"|) and 
(|87|) respectively. This will be done through a bootstrapping process and the use of 
a special energy estimate for the function w appearing in those expressions. This 
brings us to the main technical estimate of our paper which is the following: 



Proposition 5.3 (Main technical estimate). Let u = wq + w be the decomposition 
from line (192 jl . Next, let us assume that the estimate (|73|) holds with constant 2C , 
that is: 



(93) 



A 



A 2 



X(t)-2—(t)\ ^ 2Cc§—(t) + 2C(c e 
A 



A 



A 4 

sup -j(s)) A 2 . 



Then as long as the parameter X(t) is monotonically non- decreasing, one has the 
following fixed time energy type estimate for t E [0, e~ 4 ].' 



(94) 



A 



i (Arf 



1 + r 6 



(LA x w) 



(A x wf 



(t) rdr < cle 2 



A 4 

e sup -=(s) . 

0<s<t A 



Here the implicit constant depends on C from line (|93p above, but is independent of 
e and c from line (|29[) . Also, the operator A\ is defined on line (|47[) above. Lastly, 
L = dt + d r is the outgoing null derivative. 



Remark 5.4. Observe that by a direct application of estimate (|163[) of Appendix [Bl 
that the bound on line (|M|) implies: 

(95) / A- 1 -^^) rdr < c^ + e sup . 

Jm.+ l + r r o^s^t <*> 

This will be used many times in the sequel. 

We now prove the bounds (|72|) and (|73|) . This will be done separately and in re- 
verse order. To prove the second estimate ([75)1 . it will suffice for us to demonstrate 
the following set of estimates: 



Lemma 5.5 (Estimates for (173|) ). Assuming the bootstrapping estimate (|93p and 
the results of Proposition \ 5.Sl one has the following estimates where the implicit 
constant depends on line (|94[) : 



(96) \(d t u,j x )\ 3 \ < c J^+4( Coe 2 + sup ^J(s))A 2 , 

A o^ssjt A' 

(97) («, J A )A 3 = rnWH^Wy , 

(98) \(Af(u),A)\ 3 \ < e 2 ^- + ( C y + e sup ^(s))\ 2 , 
where we also have the bound^l: 

(99) lihl < e , |%| < e . 



14 The implicit constants in J99H depend only on the orbital stability bound (160 I t and thus are 
independent of the bootstrap constant C. 
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Proof of estimate (|96[) . This is the most involved of the above three estimates. To 
prove this, we begin by isolating the explicit piece involving wq from the expansion 
(f92|) . Thus, our first task is to prove that: 



(100) \(d t w J x )X 3 



< - 2 



2 



e X 



To do this we will employ the bootstrapping assumption (|93|) . We will also use the 
abstract function notation replacements from lines (|58| . Doing this, we see from 
lines (|58ap and (|58b|) that we may write: 

A A 2 

(101) (d t w ,j x )X 3 = X 4 X- 3 (FlFt) + —(FlF£) . 

Notice that from the bootstrapping assumption (|93|) and the estimate (l62j) . as well 
as the monotonicity of A we have that: 

|A| < X 2 X- X + 2C(c e 2 + sup A 4 A~ 7 (s))A 2 , 

< A^+e^ 3 < e 2 A 3 . 

Therefore, plugging this into line (|101[) and using that \(F 2 ,F X )\ < A -2 , we have 
the bound (using again (|62|)): 

\(d t w 0l jx)X 3 \ < X 4 X~ 5 + e 2 X 2 X~ 1 , 
< e 2 X 2 X~ 1 . 

This establishes ()100j) and therefore (|96|) for the wo portion of u. 
We shall now prove that: 

(102) Ud t w,J x )X 3 \ < cl^-+4(c a e 2 + sup ^(s))X 2 . 

The complication in this estimate stems from the fact that the energy estimate (|94|) 
does not provide control of the time derivative of w. In addition, (|102p is a fixed 
time estimate and therefore it is not amenable to the procedure of integrating out 
the time derivative as was done on the line ([88]) above. However we will be able to 
exploit a structure of the inner product in this expression and convert the d t w into 
LA\w derivative, which appears in J94|) . First of all, we write this as: 

(d t w,j\) = XX^ 1 {d t w,rd r Jx) . 

We now employ the following identity: 

A\{rJ x ) = 2J\ + 2rd r J\ + rA\J\ , 

along with the conditions A\J\ — and (w, J\) = to write: 

(d t w,rd r J\) = ~(d t w,A* x (rJx)) - (d t w, J x ) , 

= lidtA^w^r^^X- 1 + l([A x ,d t ]w,(rJ) x )X- 1 + XX' 1 (w,rd r J x ) . 
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Therefore, to show the estimate (|102[) we will establish the following bounds: 

(103) \(d t A x w,(rJ) x )XX 

(104) \{[A x ,dt]w,(rJ)x}\\ 

(105) \(w,rd r J x )X 2 X 

We prove these three estimates in order. By the triangle inequality we have that: 

(106) \{d t A x w,(rJ) x }\ < \{LAxw,(rJ)x)\ + \{d r Axw,(rJ)x)\ ■ 

To estimate the left hand side of (|103p involving the first term in the last sum we 
write: 

\(LAxw,(rJ) x )XX\ , 

^XXUX^Hl + rtyHAxwh^ni + r^Hr'-^F^xh^ar) , 
< XX-? \\\-?(\r)i(l + r 6 )-?LAxw\\ LHrdr) X , 

(107) < c|a 2 A" 1 + 4(c e 2 + sup A 4 A~ 7 )A 2 . 

In the last line above we have used the assumption e ^ c 2 ,. 

To conclude the estimate f|103[) . it remains to bound the second term on the 
right hand side of (|106|) . To do this we see that a simple calculation involving the 
notation on lines ([38]) and (I58b|l above, allows us to write (d r + r~ 1 )(rJ)x = Ai* 1 ^. 
This leads us to the estimate: 

\(d r Axw,(rJ)x)XX\ = |(^ A u,,F 4 )AA 2 | , 

< AAlKA^^l + r^-V^A^IU^IKl + r^^r 1 -^ 4 ^]!^^) , 

< C |A 2 A- 1 + C - | ||A^(Ar)l(l + r 5 )-ir-M A u;|| 2 2Mr) A 2 , 

(108) < cJa 2 A" 1 +cj(c e 2 + sup A 4 A~ 7 )A 2 . 

This completes our proof of the estimate (|103|) . 

To finish the proof of we need to establish the estimates (|104p - (|105p above. 
As we shall see, the proof of (|105p . with a minor exception, follows almost verbatim 
from the estimates we will use for (|104j) . Therefore we now concentrate on (|104D . 
A simple computation shows that we may write the commutator as [j4 A ,c?t] = 
— dt(Ax), which from line (|58c|) is a multiplication operator given by a function of 



s. X 1 X 
£ c^— +c,f(coe 2 + sup T7 (s))A 2 , 

A 2 A 4 

< e 2 — + (c 2 e 2 + e sup T7 (s))A 2 , 

< e ^ + ( c ^ + e sup ^( S ))A 2 . 
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the form \F%. Thus, we compute that: 

\([Ax,d t }w,(rJ)x)XX\ < \ (w, F A 12 )A 2 A| , 

< \ 2 \-^\\\-^\r)Hl + r)--ir- 2 w\\ L 2 (rdr) || (1 + r)* (r 2 ~i F 12 ) A || L2(rdr) , 

< A 2 A"I ■ (c 2 e 2 + e sup A 4 A- 7 (s))' , 

< A 4 A -5 + (cq£ 2 + e sup A 4 A- 7 (s))A 2 , 

< e 2 \ 2 \- 1 + (c?e 2 + e sup A 4 A- 7 (s))A 2 , 

where we used that AA~ 2 < e and (JTJ5]). 

The proof of the bound (|105[) is very similar to what was done above. To set 
things up in terms of the previous steps, we simply use the notation on line (|58bp 
and the Cauchy-Schwartz inequality to write: 

\(w,rd r Ja)A 2 A| 

< \ 2 \-i\\\-H\r)Hl + r)-?r- 2 w\\ L2(rdr) || (1 + r)* (r 2 -lF 4 ) A \\ LHrdr) . 
The steps are now identical to what was done in the previous due to the bound: 
(109) \\{l + r)Hr 2 -^F%\\ L , (rdr) < A -1 . 

This concludes our proof of the estimate . □ 



Proof of the identity ([97]) and the bounds ([55)1 . This follows from the orbital sta- 
bility bound (|60p and the Cauchy-Schwartz inequality. Specifically, a short calcu- 
lation shows that: 

(u,J x )X 3 = (u,(rd r J) x )XX 2 - (u,(rd r J) x )X 2 X+(u,((rd r ) 2 j) > )X 2 X . 

We leave the details of application of the Cauchy-Schwartz and the estimate 

to the reader. □ 



Proof of the estimate (|98[) . By invoking the explicit formula (|85[) and the decom- 
position it suffices to show the two estimates: 

(no) <^,Ja>a 3 < , 2 ^, 

r z A 
(HI) (^,Ja)A 3 < (c 2 e 2 + e sup ^( S ))A 2 . 



The proof of (|110p is a simple and direct calculation involving the definition (|56|) 
and the estimate (|62|) . We leave the details to the reader. 
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The proof of the second estimate (jllip follows almost immediately from (|95| . 
To see this, we compute that: 

(^,J A )A 3 < A 2 \\X-i(Xr)Hl + r)-^- 2 w\\l 2{rdr) || (1 + r)* (r 2 ^F 4 ) A \\ L ~ , 
< (c 2 e 2 + e sup A 4 A~ 7 (s))A 2 . 

This completes our proof of the estimates (|98| . □ 



Having now completed our proof of the estimate our last task in the 

section is to establish the structure estimate (1721) for the function £ defined on line 
(p?Tj) . To do this it clearly suffices to add together the following set of estimates for 
the individual terms on the right hand side of (|9Tj) : 



Lemma 5.6 (Estimates for (|72p ). Assuming the bootstrapping estimate (|93[) and 
the results of Proposition \5.3l one has the following estimates where the implicit 
constant depends on line (|M|) : 



(112) 
(113) 
(114) 
(115) 
(116) 



\{w, Ja)A| < c e 2 + c£ sup -j 

O^s^t A 

i A 4 

\{w, J\)\\ < c e 2 + c 2 sup -j 
o^s^t A 

A 2 i A 4 

\(w ,(X-2—)(rd r J)x)\ < c e 2 + c$ sup -= 

v Os£s<t A 



A 

,w ■ (2w + w) . 
( ^ . sin 



i A 4 
(2I X )J X )X\ < c e 2 + c^ sup -= 

Os£s<t A 



~ i A 4 

\{Af(u),J x )\\ < c e 2 + c 2 sup -y 

Osjs^t A 



The proof of the estimates (|112[) ~ (|116p is similar to the proof of the estimates in 
Lemma I5T51 above. We will always follow the three-step strategy: 1) Distribute cor- 
rect powers of r and (1+r) inside the inner product. 2) Apply the Cauchy-Schwartz 
inequality. 3) Refer to the estimates (JHU), & and ||73J). We will be a 

bit more terse here than before, and leave some of the details to the reader. Each 
proof will be written out under an individual heading. 



Proof of the estimates (I112p . We start with the estimate (|112l) . Using our abstract 
notation from line (|58bp we have that: 

\(w,j\)X\ , 

< A 2 A-i||A-^(Ar)^(l + r)-V 2 H| L 2 Mr) || (1 + r)^{r 2 -^F% \\ L 2 {rdr) , 

< X 2 X-i ■ {c 2 e 2 + e sup A 4 A- ? (s))^ , 

< c e 2 + c| sup A 4 A~ 7 (s) . 
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We remark here that this and similar estimates (below) are the source of our re- 
striction 4 ^ k on the homotopy class in Theorem 11.31 Notice that one cannot 
arrive at the desired bound by simply applying the orbital stability estimate (|60[) . 
It is crucial that we use (|95|) here, and this causes more weights to be placed on 
rd r J\. It is likely that one can lower the value of k in these arguments through a 
more careful analysis. We will not pursue this here. □ 



Proof of estimate f| 1 13[) . This is similar to the proof of (|112j) . Notice that from the 
estimate ([731 we have the following bound: 

|A| < e 2 A 2 + sup A 2 A- X (s) . 

Therefore, using the notation form line (|58bj) we have the chain of inequalities: 
\(w,Jx)X\ < e 2 \(w,F£)\ 2 \ + sup \ 2 \- l (s)- \(w,F$)\ , 

< e 3 + sup A 2 A^( S ) . ( c 2 e 2 + e sup A 4 A~ 7 (s))^ , 

< c e 2 + c| sup A 4 A" 7 (s) . 

Notice that we have again used the condition e ^ c§ on this last line. □ 



Proof of estimate f| 1 14[) . This will follow by a direct application of the estimate 
(175)) and the definition (f5rj|) . Notice that (fT3"|) and (|(52"|) taken together imply that: 

|A-2A a A -1 | < e 2 A 2 +c| sup A a A _1 («) . 
Therefore, a simple computation again using l|62p shows that: 

\(wq, (A — 2A 2 A~ 1 )(r<9 r J)a)| < A 2 A~ 6 • (e 2 A 2 + c| sup A a A _1 (»)) , 

< e 4 + 4 sup A 4 A~ 7 (s) , 
which is enough to imply (|114[) since e ^ c 2 ,. □ 



Proof of estimate (|115[) . The left hand side of this estimate can be bounded by the 
inequality: 

i r 2 + y- H ,^)A| < ^H^Hl + r^r-Mh irdr) 

+ || A~*(Ar)*(l + r)-2 r -2 Wo || La(rdr) . || A - *(Ar)*(l + r)-^- 2 w || L 2 Mr) . 

The estimate (|115p now follows directly from (|55|) applied to the terms involving 
u> in this last line above, and the following bound which is a consequence of the 
explicit identity ([56]) (or the notation on line l|58a[) ): 

(117) \\\-i(\r)Hl + r)-^r- 2 w \\ L 2 {rdr) < A 2 A^ . 

□ 
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Proof of the inequality (|116p . To do this, we first note that by the inequality (|86j) 
and the orbital stability estimate (|60|). used in conjunction with the Poincare type 
estimate (jSB"]) as well as the notation from line (|58b[) , we have that: 

\(Sf(u),J x )\\ < e ||A^(Ar)i(l + r)-lr- 2 u||| 2Mr) . 

By adding together the estimate (|117p and the corresponding bound (155)) for w in 
the decomposition u = wq + W, we have the single estimate for u: 

\\\-^\r)Hl + r)-ir^u\\l 2(rdr) < c 2 e 2 + sup A 4 A~ 7 ( S ) . 

Substituting this into the right hand side of the previous line we obtain the desired 
bound (|116[) . This completes our proof of Lemma 15.61 □ 



6. Space-Time Bounds and the Proof of the Main Estimate 

In this final section of the paper we prove our main technical estimate (|94[) . The 
crucial role in this will be played by the remarkable factorization property of the 
linearized Hamiltonian (|46|) . which allows us to introduce the "conjugate" Hamil- 
tonian (|50| . This new Hamiltonian H\ possesses the striking properties (f5Tj) -(|53 j) 
which are ultimately responsible for very strong estimates, proved dynamically and 
by means of simple yet quite precise physical space methods, for the corresponding 
Cauchy problem (df + H\)w — F . The key is the physical-space repulsive prop- 
erties of the operator (|50[) which lead to the desired estimates independent of how 
violently the scaling parameter A grows, so long as this growth is monotonic. This 
stands in stark contrast to the usual procedure in asymptotic stability analysis, 
which attempts to estimate the linearized operator through non- dynamic spectral 
analysis (see e.g. [S], [H]). Such a procedure is not as natural in the present context, 
which represents a truly non-linear situation not directly amenable to the standard 
perturbative techniques. From this point of view, the analysis we present here is 
close in spirit to the work of Merle- Raphael on the blow-up for the critical NLS [28] . 

We again remind the reader that it is the precise form of the non-linear equation 
(fPTjf . embodied by the first order Bogomol'nyi equation ([25]) . that is the indispens- 
able structure. 



The first thing we will need here for the proof of (|94|) is a space-time estimate 
for general solutions to the conjugated linearized equation ([50]) . For us this will 
take the form of a weighted L 2 inequality involving integration over both space and 
time variables. These are commonly referred to as Morawetz estimates, and they 
have a rich history in both linear and nonlinear analysis of the dispersive properties 
of wave equations. The estimate we use here is the based on the following energy, 
defined for sufficiently smooth and well decaying functions ip on M + : 



(118) EiM(t ,ti) 



sup 



A 



n 


A" 1 


J to JM 


+ 



(Ar) 
(1 +r 5 ) 2 r 



2 1 



(s) rdr 



(Ar) 5 V 2 
1 + r s r 3 



(s) rdr ds , 
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where < 6 <C 1 is a small fixed constant which will measure a loss in certain time 
integrations which appear in the sequel. Here we have set L = dt + d r . The main 
estimate we will use is contained in the following: 

Proposition 6.1 (Morawetz estimate for H\). Consider the time dependent Hamil- 
tonian H\. Let ip be a smooth function on [to, t±] x (0, oo), satisfying the following 
uniform bounds: 

(119) M ^ C m r , \d t iP\ + \d r i>\ ^ C m , ^ r ^ 1 . 
while decaying sufficiently rapidly at r = oo. Furthermore, suppose that: 

(120) d*ip + H x ip = d t G + H. 

Then if one has the pointwise inequalities ^5 A and AA~ 2 < e for all times to ^ 
s ^ t\, one also has the following estimate: 

(121) 

E«[V](*o,ti) < [ [ \-\Xrf [{d r Gf + e 2 {XG) 2 +H 2 ]{s) r 2 drds 

L Jt JR+ 

l + r s 

which holds with an implicit constant independent of A and 5 



sup / A ±-LjG J {8) rdr + E 5 [^](to,*o) 



Remark 6.2. The constant < 6 <C 1 will signify a small loss in time when we 
attempt to apply (|121|) in the proof of (|94|) . This is ultimately why we are restricted 
to the time interval [0, e~ 4 ] in the statement of (15. 3|) and hence in Proposition (15. 1|) . 
We also remark here that this small loss in time can in fact be avoided through 
a somewhat more careful analysis involving the precise form of the equation for A 
given on line (|69|) . 



Proof of the estimate (|12ip . Notice that all of the weights in the inequality are 
time translation invariant. Therefore, we may normalize the discussion to to = 0. 
We begin by conjugating the equation (|120p by r?. Therefore, we denote the new 
variable: 

ip = r^tp , 

We note here that the decay ip at the origin (~ r'i) will be sufficient to perform the 
integration by parts to follow. We also observe that d t ip and d r ip may be assumed 
to be bounded at r = 0. 

Next, recall that the original (super-symmetric conjugate) Hamiltonian has the 
form: 

= -d 2 --d r + V x (r) . 
r 

We define the one dimensional Hamiltonian: 

U x = -d 2 r -^ + V x {r) . 
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Then a quick computation shows that equation () 1 20(1 becomes: 
(122) d^+Hx^P = r^dtG + H) . 

The multiplier we use is the following: 



X 



l + r d 



-tL 



\~ 1+s 1 



1 



Multiplying the equation (|122j) by the quantity Xip and integrating the resulting 
expression over the interval [0,i] x (0,oo) we have the identity: 



(123) 



A 



-i (Ar 



,<5 r 



JR+ 



l + r 
fi 



L{L4,f + (Vx - - A r- 2 Wf 



dr dt 



x -iQ^ftj, G \ + d ( G) + H) -L(S) ridrdt . 
1 + r° y ' 



10 JR+ 

Here L — dt — d r is the incoming null derivative. We integrate by parts on the 
left hand side of this last expression, using the following lower bounds for terms 
involving the potential: 

c % > «-j-- 2 ) > 



.! {Xr) s k 2 
l + r 5 r 3 



which follow from (|?T|) - ([53^1 , the condition 4 ^ k, and the positivity of A. Applying 
a couple of Cauchy-Schwartz inequalities to the last two terms on the right hand 
side of (|123p . and using the positivity condition A ^ 0, we then we arrive at the 
bound: 



(124) 



i (Ar)' 



1 + r d 



^ 2 



(i) dr 



A- 1 



(Ar) a 



o Jr+ (1 + r 5 ) 2 r 



(LV) 2 drds 



l + r s r 3 



< 



A" 1 (Ar)' 5 ((9 r G) 2 +H 2 ) r 2 drdt 



(Ar) 



(1 + r s ) 2 r 



{Lxjjf dr ds 



X 



-i (Ar) 4 



R+ 



l + r' 



(m 2 + k 2 ^ 



(0) dr 



X 



i (Arf 



JM+ 



1 



L{G)-L(iP)r^drds . 



It remains to deal with the last integral on the right hand side of the above expres- 
sion. To do this, we integrate by parts with respect to the incoming derivative L. 
Employing the pointwise bound: 



< 



A 



A 2 1 + 
5, e- — —Fr? 



(Xr) s i x-1 (Ar) 5 _i 
+ X r — '—r ~ 



l + r 6 
(1 + r s )r 
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and using the equation (|122|) together with the upper bound |Va| < r 2 we have 
that: 



(125) 



(Xry 

a Jm+ 1 + r 
t 



< 



s L(G) ■ L(*fj) r*drds 
G 2 



\~ L (\r)° 



JR+ 



JR+ 



(s) r 2 dr ds 



rr7.,2 , (Ar) 5 (V-) 2 



my 



dr ds 



(1 + r 5 ) 2 r l + r s r 3 

2io 7r+ 1 + V ' Jo V 1 + ^ 



sup 



A 



-i (Arf 



1 + r" 5 



G 2 (s) rdr 



sup / A- 1 ^^^) 2 ^) 
o^ s sc*Jk+ l + r» 



Integrating by parts one more time in the term involving 9 t (G) 2 above, again using 
the fact that A 0, and using also the following fixed time Poincare type estimate: 



(126) 



A _1 (Ar) 5 G 2 dr 



< 



\-\\r) s (d r G) 2 r 2 dr , 



we add together the estimates (|124[) - (|126[) and take the sup over different times to 
achieve the bound: 



(127) E g [i>](Q,t) < 8- 1 



JR+ 

sup 



X-'iXr) 5 [{d r Gf + e 2 (AG) 2 + H 2 ] (s) r 2 drds 
A 1 -=-^G 2 (s) rdr + I tf [#|(0,0) 



i (Ar)« 



where we define the one dimensional energy analogous to (| 1 1 8[) : 

-i (Arf 



E 5 [V](0,i) = sup / A" 

o^ssst Jk+ 1 + r a 



(LV) + fr 



(s) dr 



A - 



(AO 5 , r T, 2 , , 2 (Ar) 5 ^ 2 



1 + r" 5 r 3 



(s) dr (is . 



(1 + r 5 ) 2 r 

Finally, to complete the proof, we use the expansion: 

r^^Lip — Lij, H — r _1 ?/> 



and the fact that 4 ^ k to bound the energy Ea from above and below by E$ . This 
completes our proof of the estimate (] 1 2 1 [) . □ 



We now turn to the proof of (|94|) . The precise statement of what we need to 
show is the following: 

Proposition 6.3 (Energy estimates for the quantity w). Let u = Wq + w be the 

decomposition of u given on line (I92[) , where u itself is part of the decomposition 
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(|44| of the full field <j>. In particular u solves the equation (|45|) . Suppose also that 
the initial conditions for cj> are given as to satisfy ()28|) - ()29j) . and that u obeys the 
extra decay estimate (|147p (this will be proved shortly). Furthermore, assume that 
the main assumptions of Proposition \5.3\ hold, in particular we have that (|93[) and 
A 0. Then the following estimate holds for t € [0, e -4 ].' 

(128) E s [A x w]{0,t) < 5- l (cle 2 + e sup ^(s)) . 



The implicit constant depends on C from line (|93[) but is independent of the size of 
c 0; e, or 5. 

In particular, for a fixed < S <C 1 we have that the estimate (|94p holds. 



Proof of the estimate (|128[) . The first order of business is to reduce the proof to 
simpler bounds. In the sequel, we will only show that: 

/ A 4 \ 

(129) Ej[A A ti;](to,ii) < S- 1 (E s [A x w}(to,t ) + e 3 + e sup -=(«)) , 

for all time intervals [io,ti] inside the regular interval [0, T**), where again T** is 
the blowup time, provided that one also has the inequality: 

(130) f 1 ^J(s) < e . 

•/to A 

We claim that along with the bootstrapping assumption (|93|) and the analysis done 
in Section [O] this is enough to establish (I128j) . 

To verify this claim, first notice that if we are in the time interval [0, T*] where 
T* is defined as on line (|78|) , then we automatically have (| 1 30[) on account of (|79|) . 
Therefore, we may work inside intervals of the form [T* , t], and we are only trying 
to establish: 

(131) E s [A x w](T*,t) < S-h^it) . 

Notice that we are using the monotonicity of A 4 A -7 established in Section I5TT1 We 
now claim that (|131[) easily follows from (| 129[) and the bootstrapping assumption 
Indeed, let [tfe_i,tfe] be any interval where equality in (| 1 30[) holds. Then we 

have: 

(132) / ^-H\ 4 \- 7 )ds > ±ds > e- 3 Tvds = e~ 2 . 

Here we have used the bootstrapping bound (|93| in the simple form |A — 2A 2 A~ 4 | <C 
A 2 A _1 , which holds as long as we are in the region past [0,T*] (in particular, one 
has access to a lower bound consistent with (|78j) which allows one to uniformizc 
the RHS of (f9"3"]l ). Notice that we have also used the orbital stability bound ([52^1 
several times in deriving the inequalities. Integrating the inequality (|132[) we see 
that: 
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on any time interval [tk-i, tk] past [0, T*\ where (|130j) also holds. It is now a simple 
matter to derive (I131j) from (I129j) . We first decompose the interval [T* , t] into a 
finite collection of N subintervals [tk-i,tk] where equality in (|130[) holdsj3 On 
each of these intervals, we have the estimate (|129[) . By repeatedly using the bound 
(|133|) . each of these estimates may be inductively expanded to yield: 



k 



E s [A x w](t k ^,t k ) < J-iC^f— J -ejj(t h ) 

, 5~ l C A 4 . . 
~ <^- e A^ W ' 

where C\ ~ e e is some incredibly large constant that beats the (uniform) im- 
plicit constant C<5 -1 appearing in the estimates (|129[) . Summing this last line over 
k < N, we have the bound fl3"T]) . 

We now prove (| 1 29[) under the additional assumption that (|130[) also holds. We 
start by providing the general setup, and then reduce the proof to a number of 
separate estimates to be dealt with under their own bold-faced headings. We first 
record the equation for w. Recall that the purpose of the decomposition (|9"2"|) is 
to eliminate the main source term A X {I X ) on the right hand side of (|4"5|) , obtained 
after applying A x . Therefore, we have that: 

A x [d 2 w + H x w] = -A x (w Q ) + A x Af(u) . 

To put things in the form where the estimate (|12ip can be used, we commute the 
A\ operator with df on the left hand side of this last equation, which yields the 
expression: 

d 2 t W + H X W , 
= - A x (w ) + [<9 t 2 , A x ]w + A x Af{u) , 

- d t (A x )d t (wo) - d t (A x d t (w )) + 2d t (d t (A x ) ■ w) - d 2 t {A x ) ■ w + A x N{u) , 
(134) 

= Mi + d t M 2 + d t Ri + i? 2 + R 3 , 

where we have set W — A x w. The M terms above constitute the "main source" 
which feeds the quantity W through the wave-flow of the Hamiltonian H x . By 
contrast, the R terms on line (|134p are for the most part "errors" which will be 
reabsorbed back onto the left hand side of the estimate (|128p . This is where the 
limits on the time interval and the decay estimate (| 14T[) will come in to play. We 
now turn to the details of all of this. 



12 On the last interval there may be a strict inequality in (1 1301 1 ■ but this single interval may 
also be estimated with l|129j l. and the answer may then be directly added into the final bound. 
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To estimate the terms Mi and Mi via the general bound (|121[) . we will show the 
following four estimates: 

(135) I' f A- 1 (Ar) 4 |Mi| 2 (s) r 2 drds < e 3 + e sup ^(s) , 

Jt JR+ to^s^ti A 

41 ' A 4 



(136) / / \- l (\r) s \d r M 2 \ 2 (s) r 2 drds < e 3 + e sup -y (s) , 

(137) / A-^Ar^AMal 2 ^) r 2 drds < e 3 + e sup ^(s) , 

./to JR+ to^s^ti A 

/" fAr) 5 A 4 

(138) sup / A- 1 -^ T |M 2 | 2 (s) rdr < e 3 + e sup — (s) . 

Recall we are assuming that e ^ c 2 ,, so these estimates will be enough to generate 
the right hand side of (|128p for the M terms. 

We shall estimate the R terms on line (|134jl in a nonlinear fashion. Specifically, 
we will bound them in terms of a small constant times the energy on the left hand 
side of (|128p . plus one term involving the nonlinearity J\f applied to wo, which fits 
into the pattern of the right hand side of (|135p - (|138[) above. What we propose to 
show is the following: 

(139) 

/ A- 1 (Ar) 4 |9 r J?i| 2 (s) r 2 drds < e 2 E s [W](t , h) , 

to JM+ 

(140) 

r*i r 

/ A- 1 (Ar)*|Afli| 2 (s) r 2 drds < e 2 E s [W](t , h) , 
(141) 

sup / A- 1 -^^!! 2 ^) rdr < e 2 E 5 [W](t , *i) , 

t ^s<ti JR+ 1 + T 

(142) [ A- 1 (Ar) [5 |i? 2 | 2 ( S ) r 2 drds < e 2 E s [W](t ,h) , 



Jt JR+ 

(143) I' I A- 1 (Ar) 5 |i? 3 | 2 (s) r 2 drds < e sup ^J(s) + t{t 2 E s [W]{t^ h) . 

Jtg JR+ to^s^ti A 

Using the conditions e is sufficiently small, and that t\ ^ e~ 4 , all of the estimates 
added together will imply the estimate (|128j) for the R terms on line (I134| above. 

We now turn to the details of the proofs of the estimates ()135|) - (1138|) and ()139l) 
(|143p . We will do each of these separately and in order. 

In what follows, we will consistently use the following "identities" which are in 
accordance with our notation from Section 12.11 

where the second "inequality" holds provided that A ^ 1. 
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Proof of estimate (|135|) . Here and throughout the sequel we will rely heavily on 
the abstract function notation from lines |58|) above. Multiplying together estimates 
from lines (|58ap and (|58cp we have that: 



\d t (A x )dt(w )\ < (A 4 A- 5 + |A|A 2 A- 



F 9 X 



To resolve the second term above which contains the expression A we do the fol- 
lowing. Notice that the orbital stability bound AA~ 2 < e and the bootstrapping 
estimate (|93p give the rough pointwise bound: 

|A| < A^ + ^ + A^-^A 2 , 
(144) < A 2 A^ +e 2 A 2 . 

Therefore, substituting this estimate back into the previous line we have that: 

\d t (A x )d t (w )\ < (A 4 A- 5 + e 2 A 2 A- 2 )F| . 
Using this last line we can now estimate: 

(L.H.S.)GE& < P f A- 4 (a 4 A~ 5 + e 2 A 2 A- 2 ) 2 ^ r 2 drds , 

< f 1 A 8 A- 14 (s) ds + e 4 f 1 A 4 A- 8 (.s) ds , 

Jt Jt 

< sup A 4 A~ 7 (s)- f 1 X 4 X- 7 {s)ds + e 4 sup A 3 A~ 6 (s) • f * AA~ 2 (s) ds 

to<s^ti Jto fo^s<ii Jto 

< e sup A 4 A~ 7 (s) + e 7 . 

To obtain the last line, we have used both the estimate (|130[) . the bound (|62p . and 
the assumption that A 1. 



Proof of estimate (|136[) . This is very similar to the analysis above, with an 
addition of a small twist. First of all, by combining lines (|58a[) . (|58c[) . and 
then (|144p we have the abstract notational bound: 

\d r (A x dt(w ))\ < (a 3 a- 3 + |a|aa- 2 )f 4 , 

< (A 3 A- 3 + e 2 A)f 4 . 
Substituting this into the left hand side of (|136|) we have the chain of inequalities: 

(L.H.S.)(H3i < ^ f A- 4 (a 3 A- 3 + e 2 A) 2 F 4 r 2 drds , 



< 



41 A 6 A- 10 (s) ds + e 4 ^ X 2 X- 4 (s) ds , 

to J to 

ti rti 



< f 1 A 6 A- 10 (s) ds + e 4 sup AA- 2 (s) ■ f ' AA~ 2 (s) ds 

Jta to^s^ti Jto 

A 6 A" 10 (s) ds + e 5 . 
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We now need to do a little work, because the first term on the right hand side of 
this last line above is not manifestly of the correct form. To correct it, wc first 
integrate by parts with respect to time which yields the identity: 

f 1 \ 6 \- w ds = A 5 A- 9 (ii)-A 5 A- 9 (£ ) + 10 f 1 \ 6 \- 10 ds-5 V \\ 4 \- 9 ds . 

Jto Jto Jto 

Now, using the bootstrapping assumption (|93|) as well as the estimate f|62[) . this 
last expression leads to the following nonlinear bound: 

(145) 

41 A 6 A- 10 ds , 

<c| / 1 A 6 A~ 10 (s) ds + e sup A 4 A~ 7 (s) + f * (e 2 + sup A 4 A~ 7 ) • A 4 A~ 7 (s) ds 

Jt *o^s^*i Jt to^-^s 

< 4 f 1 A 6 A~ 10 (s) ds + e sup A 4 A~ 7 (s) + e 2 f * \ 4 \- 7 (s) ds , 

Jto to^S^t! Jt 

ftl 



\ b \- w {s)ds + e sup A 4 A- 7 (.s)+e J 

to ^S^tl 



To go from the first to the second line above, we have used ([75]) and the monotonic- 
ity established in Section [5] (this works as long as T* ^ to, whereas in the other 
case we may as well assume that to — 0). Notice also that in the last two lines 
above we made several uses of the assumption (|130l) . 



Proof of estimate (|137[) . This is virtually identical to the proof of (|136[) . Another 
simple calculation using lines (|58a|) and (|58c|) . and then (I144|) gives us: 

\X(A x d t (w ))\ < (A 3 A- 3 + |A|AA- 2 )^ A 3 , 

< (a 3 a- 3 + £ 2 a)f 3 . 

Substituting this into the left hand side of (|137[) the proof follows verbatim from 
the calculations done in the previous paragraph. 

Proof of estimate (| 138[) . Once again using lines (|58a[) and (|58c[) . and then (|144[) 
we have that: 

\(A x d t (w ))\ < (a 3 a- 4 + |a|aa- 3 )^ 3 , 

< (A 3 A- 4 + e 2 AA- 1 )^ 3 . 

Plugging this last line into the left hand side of (|138[) and simply using the bound 
(|62p we arrive at the chain of inequalities: 

(L.H.S.)OSS) < sup / A" 1 (a 3 A~ 4 + e 2 X\- 1 ) V 4 rdr , 

t ^s^t 1 Jm+ ^ ' 

< sup A 6 A~ n + e 4 sup A 2 A" 5 , 

< e 2 sup A 4 A" 7 + e 6 . 
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We now turn our attention to the proofs of the estimates (|139p ~ (1143p . 



Proof of (|139[) . We first provide a pointwise bound for the term d r R\ from line 
(|134p above. This involves a simple application of the abstract notation (|58cp , the 
rules from line (|38p . and the decomposition: 



d r w — — W H — cos(I\) ■ w . 
r 

Together, these give us the following estimate: 

\d r (dt(Ax) ■ w)\ < (\XXw\ + \XW\)F% . 



Plugging this identity in the left hand side of (| 139[) leads us to the estimates (also 
using (|62|)): 



to JR+ 
ti r \2 



< 



Jt J&+ A 5 



X~ 1 (Xr) 5 \d r Ri\ 2 r 2 drds 
;2 W 2 



((r 4 + r e )F 12 ) x rdrdt 



< f 2 



-i (Ar) 



s r w 2 



to 
ti 



(1 +r) s 



W 2 



rdr dt 



_! {\r) s W 2 



< f 2 I I X- 1 rdrdt 



<e 2 



'to JR+ (l+r) s r 3 
E s [W}(t ,h) . 



To obtain the second to last line above, we have used the comparison estimate (|164p 
from Appendix IB1 on the term involving w. 



Proof of estimate (|140p . This is virtually identical to the proof of (| 1 39[) in the 
previous paragraph. A simple calculation using the notation from lines (|58[1 and 
line ([38]) gives us the bound: 



\X(d t (A x ) ■ w)\ < XX\w\F 7 x . 
The proof now follows line for line from the calculations performed above. 

Proof of estimate (|14ip . We again use the formulas on line (|58cp which give us: 



\(dt(A x ) 



< X\w\ F{ 
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Substituting this last line in the left hand side of (|14ip . we have the following chain 
of inequalities where the second to last line involves the bound (|164l) : 

sup / X -^m»r*r< sup / £ (r^) x rdr dt , 



> _ 1 (Xr) 5 w 2 



^ 6 sup / A — — — rdr at 



< e 2 sup / A t~ — — —rdrdt 



< e 2 Etf[W^(t 0) ti) ■ 



(l + r)' 5 r 2 



Proof of estimate (|142p . We estimate the i?2 term from the line (|134|) . By using 
the second derivative identity from line (I58c|) . the estimate (|144p . and then the 
bound (|62p we have the following pointwise estimate: 

\d 2 (A x )-w\ < (X 2 \-'+ ( 2 \ 2 )Fl-\w\ , 

< e 2 \*Fl ■ \w\ . 

We now substitute this estimate into the left hand side of (|142p which allows us to 
estimate: 

/ 1 / A -1 (Ar) ff |i?2| a r 2 drds < e 4 f'f X 5 {Xr) s \w\ 2 ■ r 2 drds , 
Jto JK+ Jto JR+ 

< e 4 f 1 [ X-HXrY^C -{r G F u ) x rdrds , 

Jto Jr+ r 

41 f X-i-^^rdrds 
41 ' _! {Xr) s W 2 



X — -j- — — rdrds , 

/to JR+ (± + r) s r3 

< e 4 E 5 [W](* ,ti) • 



Proof of estimate (|143p . First of all, using the formula for the nonlinearity Af(u) 
given on line (|85p . and by making use of the formula (|4"T|) for the operator A \ as 
well as the formula (|56p for wq, we easily have the pointwise bound: 

(146) |^(.)| < x^Fl + (H + y-H + tN±M . 

We will deal with the first term on the right hand side above by itself. The other 
two terms can be handled together. 
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We now substitute the first term on the right hand side of (I146|) for R3 on the 
left hand side of (|143|) . Doing this we are left with estimating: 



in 



\"\- lL {\r) s Fl 2 r 2 drds 



< 



A M A- i4 (s) ds 



< sup A 4 A~ 7 (s) 

< e sup A 4 A" 7 (s) 

io ^s^ti 

This proves the estimate (|143|) for the wq portion of R3 . 



\ 4 \- 7 {s) ds , 



in 



It remains to deal with (I143|) for the last two terms on the right hand side of 
(|146j) . Upon substitution of these into the right hand side of (|143j) we have that: 



< 



< 



t] 



A-^Arf (M 2 + K| 2 ) 
,-1 (Ar) 



r-4 



r 2 dr ds 



\S r w 2 



sup (1 + r) d \u\ 2 + e 2 

r 



W/2 



(1 + r) 5 (\u\ 2 + \w \ 2 ) rdrds , 



-1 (Ar) 



S r w 2 



to JR+ 



(1 + rf 



W 2 



rdr ds , 



< tU 2 



tl 



X-^-^LK rdrds 



i+ (1 + r) s r 3 
< 4e 2 E s [W](to,h) . 

Notice that in the above estimates we have made crucial use of the special pointwise 
estimate (|147l) proved below. This is the only place in the paper which requires the 
extra decay of the initial data. This completes our proof of the estimate (|143[) . and 
thus our proof of Proposition 16. 31 □ 



6.1. A Simple Decay Estimate. In this subsection we will prove the rough decay 
estimate: 

(147) sup {l + r) 5 \u\ 2 < t s e 2 , 

r 

That is, our aim is to show that the reduced field quantity u enjoys some amount 
of pointwise decay outside of a sufficiently large cone centered at the space-time 
origin t — and r = 0. 



Lemma 6.4 (Decay of u at space- like infinity). Let u = <f) — I\ be the reduced field 
quantity as defined in Lemma \4-.1\ which in addition satisfies the initial conditions 
of Theorem \l-4\ In particular, u is a solution to the equation (|45| with initial data 
(f28|) -(|29 |) and obeys the estimate ([60]) on the time interval [0, T] where <p exists and 
remains smooth. Then u also obeys the following stronger energy type estimate for 
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any time t € [0,T], for which in addition A 1: 
(148) / 

Remark 6.5. To transform estimate ()148|) into an L°° bound can be done in an ele- 
mentary way by applying the Poincare type estimate ()68|) to the quantity r\3t^rU, 
and then using the bound (I148j) to estimate the resulting right hand side. Here 
X3t^r is a smooth cutoff onto the region where 3t ^ r which satisfies the homogene- 
ity bound |X3{<rl ~ r~ x . Therefore, we arrive at the estimate: sup 3t ^ r r|u| < e. 
By combining this with the pointwise bound \u\ ^ e which holds everywhere, we 
easily have (|147p whenever 6 ^ 1. 



(dtc/l) 2 + (d rU ) 2 + XrU 



rdr < e . 



Proof of estimate (| 148[) . The proof is an integration by parts argument with a cer- 
tain multiplier. We denote by ct{y) a smooth increasing function, supported where 
10 < y, satisfying a' ^ 3y and the homogeneity bound y~ x a ^ a'. The desired 
result will now follow from computing the left hand side of the identity: 



(149) 



d 2 cj) + H x u- Af(u) d t <p ■ a(r - 2s) rdr ds = 0, 



where Af(u) = R.H.S. (|4"5"]) . Also, we will write the Hamiltonian from line (|46|) as 
H\ = —d 2 — r~ 1 d r + Q\. Notice that we have Qx ^ ck 2 r~ 2 on the support of 
a(r — 2s). 



Using the factorization (|4"6"| as well as the fact that Ax(Ix) = 0, thanks to 
we may transform (1149)) into the identity: 



(150) 



(C^) 2 + {d r uf + Q X U 2 



rdr 



(d t <j>) 2 + d t <j)d r u + (d r u) 2 + Q x u 2 ■ a'(r - 2s) rdr ds 



cos(Ix) Ix • o! rdr ds — — 
o Jm+ r 2 
t 

J\f(u) d t 4> ■ a rdr ds , 



Qx u ■ a rdr ds 



= T 1 +T 2 +T 3 +T 4 . 

The proof will be complete once we show that the terms on the right hand side 
of this last expression are either non-negative or are bounded in absolute value by 
Ce 2 . In fact, it is more or less immediate that we have: 



(151) 



< 



r!+T3+T 4 > 



The first estimate above is a consequence of the Cauchy-Schwartz inequality, the 
orbital stability bound (|60|) . and the following fixed time estimate valid for 1 ^ A: 



\h ■ ct || i 2( rdr) 



< 
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This last line uses our assumptions that 4 ^ k and 1 ^ A. Specifically, the 
ODE bound (JHU) | A A 1 1 < eA (from (JHU) above) and a simple calculation, us- 
ing the assumption that 1 A and involving lines (|48[) and (|34p, give us the bound 
\h-a'\ < \ri\ <e(l + r)- 3 . 

The second bound on line (|151| will follow from the estimate \T^ + T±\ < eT±. 
The desired result is then a consequence of the homogeneity property of a and 
bounds: 

Mu)\ < eM , \Q X \ < el . 

The first bound above is a simple consequence of the orbital stability estimate ([SO]) 
together with (ft?5|) . The second bound follows again from the estimate | A A 1 1 < eA 
of (jnij) and the explicit formulas on lines (|34p - (|3"S"|) . 

The estimate (I148[) now follows from the form of the left hand side of (|150|) and 
the smallness condition (|bu]) . □ 



Appendix A. Computation of the constant C* 

The purpose of the appendix is to derive an explicit formula C* = for the 
special constant C* which appeared on line ([50)) . Here we have written J = J\ 
according to previous notation. In what follows we shall also denote I = I\ = I k . 
Rescaling we have that: 

C» = Ti + T 2 + T 3 , 

where: 

9 /(aJ + br 2 j) 2 , \ 
Ti = - k 2 ^ >- , sin(27) • j) , 

T 2 = (aJ + br 2 J , rd r /> , 

T 3 = - (rd r (aJ + br 2 j) , rd r j\ . 

Recall that the constants a and b are given on line (|57p . Using now the identity 
sin(2J) = 2 sin(7) cos(J) as well as (|48|) and (|25p . we have that: 

T 2 = -k 2 I (aJ+ br 2 J) 2 d r ( sin 2 (J)) dr , 

= - / (a + 6r 2 ) 2 J 2 9,.(J 2 ) dr , 

(152) = 2ab / J 4 rrfr + 26 2 / J 4 r 3 dr , 

To obtain the last line above, we have used the expansion J 4 — fc 2 J 2 sin 2 (J), the 
Pythagorean identity, and the definitions of a, b. 
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We now move on the term Ti above. Here we have directly that: 
T 2 = f (aJ + br 2 J)d r J r 2 dr , 

JR+ 

= -a I J 2 rdr -2b f J 2 r 3 dr , 
Jr+ Jr+ 

(153) = - 7 / J 2 r 3 dr ■ 

4 Jr+ 

Finally, we compute that: 

T 3 = - a I {rd r J) 2 rdr - b j (rd r J) 2 r 3 dr - 2b I Jd r J r 4 dr , 
Jr+ Jr+ Jr+ 

= - ak 2 j J 2 cos 2 (7) rdr - bk 2 [ J 2 cos 2 (7) r 3 dr + 4b [ J 2 r 3 dr 

JR+ JR+ JR+ 

(154) 



a I J 4 rdr + b / J 4 r 3 dr + 46 / J 2 r 3 dr . 

JR+ 



We now add together lines (]152|) (|154[) into the single formula: 

(155) C* = \a ( J 4 rdr + \b [ J 4 r 3 dr + 36 f J 2 r 3 dr . 

2 Jr+ 2 J R + J R + 

It remains to compute the first two integrals in this last expression. 

J+rdr = -kf rd r (cos(I))j 2 rdr, 
Jr+ 

= 2k j cos(/)J 2 rdr + 2k 2 [ cos 2 (/)J 2 rdr , 
Jr+ Jr+ 

= k 2 / r5 r (sin 2 (/)) rdr + 2k 2 / J 2 rdr - 2 J 4 rdr , 

Jr+ Jr+ Jr+ 

= (2k 2 -2) I J 2 rdr -2 / J 4 rdr . 

JR+ JR+ 

An almost identical calculation also shows that: 

J 4 r 3 dr = {2k 2 - 5) / J 2 r 3 dr - 2 / J 4 r 3 dr . 

Therefore, recalling the definition of a and 6, these last two calculations together 
give: 

a f J 4 rdr + b [ J 4 r 3 dr = --b [ J 2 r 3 dr . 
Jr+ Jm+ 3 J R + 

Inserting the last line into (| 1 55|) and using that 6 = i we have C* = as desired. 
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Appendix B. A general functional analysis lemma 

In this appendix, we prove a general form of a coercive estimate we need through- 
out the paper. This turns out to be more expedient, because the required struc- 
ture is simply a matter of compactness and weak convergence in various weighted 
Sobolev spaces. The general result which we propose to prove here is the following: 

Lemma B.l (Coercive bounds for first order operators). Let Bi be a sequence 
of first order differential operators with real smooth (but not necessarily bounded!) 
coefficients on the half line (0, oo), continuously indexed (in the weighted L 2 space 
defined by the LHS of (|156|1 below) by £ £ [0, 1] and such that the following sub- 
coercivity holds for some continuously indexed (for functions in the norm (|157[) ) 
function \f t \ < r - 2 -T+ CT (l + r yi° w uh a > 0: 

(156) J r+ rdr = | + [_ + + M ] rdr , 

for any real valued function ip with finite norm: 

Here ^ 7 is a fixed parameter, and < C ^ hi is some strictly positive function. 
Then there exists a universal constant, uniform in £, such that the following bound 
holds: 

(158) < I ( -^rdr, 
for any real H, 1 function ip which also satisfies: 

(159) / V • J 1 m e rdr = , 

Jr+ 

for some positive weight function < mg. Here the function J 1 is the (nontrivial) 
"ground-state" given by BgJ e — 0, and we are assuming mi J £ (7i 7 )* depends 
continuously on £. 

Proof of Lemma \B.1[ The proof is based on a contradiction argument centered 
around weak convergence. Suppose that the estimate (| 1 58[) was not true. Then 
there would exist a sequence of tp n and £ n such that: 

(160) f {Bf J ' n)2 rdr < en ||^||^ , 

where c„ — * is some sequence of constants. We assume that this sequence is 
normalized so that || ip n ||^ 7 = 1. The space TL 1 is a Hilbert space (with an ob- 
vious scalar product) defined as a closure of Cq°(R+) functions in the Ji 1 norm. 
Therefore, we can choose a subsequence ip nk which converges weakly in H, 1 to 
?/>oo € TT' . Furthermore, we may assume (by perhaps taking another subsequence) 
that £ nk — > £oo for some £oo G [0, 1]. We now use ifj n and £ n to denote this subse- 
quence. Also, note that by Cauchy-Schwartz the unit normalization implies that 

ll^oollw < I- 
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By the continuity of the J 1 , and the uniform boundedness of the ip n , we have 
from the identity (I159P that the limiting function satisfies: 

(161) / Voo • J £o ° m ioo rdr = . 

Jr+ 

Therefore, since the "ground-state" J e °° is unique (it satisfies a first order ODE) 
and the measure mg^rdr is strictly positive on (0, oo), we will have a contradiction 
if we can establish that i/'oo is nontrivial. This contradiction would come from again 
invoking uniform boundedness and the continuity of the operators Bg which implies 
that: 



rdr — ► 



so that ipoo is a weak solution of Bg^ipoo — 0, and hence a smooth solution via ODE 
regularity, thus violating uniqueness as (|16ip implies ipoo =^ f3J e ' x for any constant 

To show that ip^ is nontrivial, we make crucial use of the sub-coercivity condi- 
tion (|156p . By the unit normalization, the universal lower bound on hi, and the 
assumption that (|160j) . we have that there exists a universal lower bound to the 
limit: 

> lim / fe n ipn rdr . 



n — >oo 



+ 



Therefore we shall have that V'oo is not everywhere zero if we can show that the 
sequence fi„ipn converges strongly in ^(rdr). This in turn follows from the uni- 
versal bounds on and continuity of fg, and fact that i\) n converges strongly in the 
weighted space: 

r- 2 -T +CT (l + r)- 2CT ^ 2 rdr = \\ ^ \\% ^ . 

This latter strong convergence is provided via uniform boundedness and the com- 
pact inclusion TL 1 (s £ 2 CT whenever < a. □ 

In practice, we will only need two special cases of the Lemma IB . 1 1 above. The 
first case is where Bg = A\ and 7 = 0, where A\ is the first order operator from 
line (l47l) above. The second cases are when we set £ = A -1 , with 1 ^ A and: 



. » 



B t = (1 + (A-V)) 5 A, (1 + (A-V)) 5 , 

where in this case we set tjj = (1 + (A _1 r)) _ ^u, as well as J e = (l + (A _1 r)) 2 J\. 

Finally, in this case we set me = (l + (A^V)) 5 . We apply this to 7 = 2 — S and 
7 = 3-5. 

In all of these cases we leave it to the reader to prove that the condition (|156|) 
holds (the continuity is obvious). This is a simple matter of integration by parts, 
the fact that 4 ^ k (notice that this works for our range of 7, which is the main 
thing to check here), and also that we have chosen 5 <C 1. However, we do call 
the readers attention to an important and perhaps subtle point. In order for the 
integration by parts to work, it is necessary to show that a boundary term of the 
form lim r ^o r _7 -0 2 vanishes for any function ip in the space H? ' . This follows from 
the finiteness of that norm, and the fact that the Poincare type estimate (|68)) above 
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applied to r 2 ip implies that this function is continuous on the closed interval [0,1]. 
This latter fact is perhaps a bit subtle, and it is crucial for showing r~irb vanishes 
at r = via the finiteness of the weighted L 2 norm (no derivative) contained in 
TL 1 . Again, we leave the reader to check the details of all this. 

Now, Applying the above result in these two cases and rescaling by A, we have 
that: 



Lemma B.2 (Applied version of Lemma |B.1[) . Suppose that 4 ^ k and A\ is the 

operator defined on line (|47[) . Then ifu is a function which satisfies the admissibility 
condition of Lemma \1.9l and is such that the energy from line (160[) is finite, and 
such that the orthogonality relation holds: 

u ■ J\ rdr = , 

then one had the following universal bounds whenever 5 <C 1 is small enough ( and 
1 $C A in the last two cases): 



rdr < I rdr 

{l + r) s r 2 



(162) 


JI+ r 


(163) 


f (Xr) s u 2 
J R+ (l + r)*r4 


(164) 


f (Xr) 5 u 2 
V {l + r) s r 5 



+ 



(i + r y r 3 
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